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ABSTBiCT 


^ A simulation programme is developed to predict the profiles 
of liquor temperature, concentration, heat flux, overall heat transfer 
co-efficient, pressure, pressure gradient, void fraction and vapor quality 
along the tube length of a natural circulation LTV evaporator. These 
quantities are predicted for the nonboiling and boiling sections of an 
evaporator tube from a ccmbination of the following expressions scontinuity, 
momentum and energy balance equations, and frictional pressure drop and 
holdup correlations for^ sin^e plase and two phase flow, heat transfer 
ooirrelation for the nonboiling, nucleate boiling and annular film regions, i 

and correlations for the engineering properties of the process liquor such : 
as thermal conductivity, viscosity, density, surface tension, boiling point | 

rise and specific heat, A stepwise computation procedure has been developed ! 

starting with feed inlet conditions assuming a series of differential sections 
in the nonboiling and boiling sections of the tube and using the criteria of 
a maximum in liquor temperature profile and a minima in the pressure gradient 
curve for the transition from nonboiling to nucleate boiling and from nucleate 
boiling to annular film region respectively. 

Results from this simulation progcamme agree with Badger’s 
experimental data with an overall accuracy of 15^ and agreement is within 
396 for the temperature profile and an average deviation of 25?^ for the 
overaj(.l heat transfer co-efficient based on Gudmundson's experimaital data 

*m^ ■ ■ ^ f 

for water » The model is applied to bamboo, baggise and pine black liquors i \ 
over the concentration range 12-28^ and the pronounced effect of viscosity 



on heat transfer rates is illustrated for the f-oUowing conditions s 

feed-rate 100 tons/hr, temperature 50° C and concentration 12.0^; 

Evaporator-tube diameter 51 mm, length 7*6 m and number of tubes 500, 

dirt scale resistance = 0*45 m °K/KW. The nonboiling section length 

is 0*92 m and 2.15 m for bamboo and bagasse black liquors; the corresponding 

values of the overall heat transfer coefficient are 2. 4-2. 7 KW/m “K 

(from 50°C, 0.9 cp to 106°C, 0.45 cp) and 0.8-1^1 KW/m^ °K (from 50°C, 

4.7 cp to 102. 4°G, 1.9 cp) with clean surfaces. In the boiling section 

the overall heat transfer co-efficient values are 2.4-5»1 "K and 

2.5-5 *5 KW/m °K for bamboo and baggase black liquors respectively; the 

final concentration, tempera tvire and viscosity values are 95.5°C, 

1.1 cp for bamboo and 25‘9?^» 86.2°C, 9«8 cp for ba^sse black liquors 

assuming clean tube surface. With, the inclusion of the dirt scale 

resistance the nonboiling length increases to 1.85 m and 2.75 m with a 

'2 

corresponding decrease in heat transfer 00 -efficient to 1.2-1 .5 KW /m °K 
and 0.6-0 .7 KW/m °K for bamboo and baggase black liquors respectively in 
the nonboiling section. In the boiling region the dirt resistance 
decreases heat transfer co-efficient to 1.2-1 .4 KW/m °K and 1 .2-1 .45 KW/iii ' 
for bamboo and baggase black liquors respectively giving a final liquor 
with ccxicentration, temperature, viscosity of 15 105° C* 0*41 cp for 
bamboo and 15*49^» 99°C, 2.0 cp for baggase black liquor. 
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CmPTER 1 
INTRODUCTION 

Natural circulation long tube "vertical (LTV) e'vaporators are 
widely used for the concentration of solutions such as milk, fruit 
juice, spent pulping liquors etc. This study deals with the develc®’’* 
ment of a model to predict changes in e'vaporation capacity and to 
recommend auxiliary heat exchange equipment that may be needed to acco- 
modate a large change in the viscosity of the process liquor. A good 
example is the e'vaporation of weak black liquors from alkaline pulping 
process. in.per mills in India use bamboo as the major raw material and 
smaller amounts of mixed hardwoods, salai, bag&se, grasses, pine etc. 
Increasing quantities of bag^se and eucalyptus may also become ffivail- 
able to the paper industry. The evaporation plant in a p'aP®^ mill must 
be capable of handling bamboo, bag^se, eucalyptus black liquors and 
also mixed black liquors whose viscosities "vary considerably over the 
same range of concentration and temperaturi^^^ . • 

Fluid flow and heat transfer mechanisms in a LTV e"vaporator result, 
in nonboiling, nucleate boiling and annular film regions along the 
tube length and determine the changes in solution tsnperature and 
concentration. In simulating the performance of a LTV evaporator process 
feed liquor conditions such as flow rate, temperature, concentration and 
pressure are usually available. Simulation consists in predicting the 
profiles of liquor temperature, concentration, heat flux and overall 
heat transfer co-efficient along the length of the evaporation tube. 

These computations would also give variation in pressure, pressure 
gradient, void fraction and "vapor quali'ty along the tube length. All these 
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quantities can be predicted for the given e^porator unit from a combinatic^ 
of the following expressions; Continuity, momentum and energy balance 
equations, and frictional pressure drop and holdup correlations for 
single phase and two phase flow, heat transfer correlation fcr the 
nonboiling, nucleate boiling and annular film regions, and correlations 
for the engineering properties of the process liquor such as thermal 
conductivity, viscosity, density, surface tension boiling point ^lise 
and specific heat. Computations would begin with the feed inlet condi- 
tions at the bottom of the evaporation tube and proceed through the nonb- 
oiling, nucleate boiling and annular film regions, and thus tmverse 
the entire tube length. Computations are carriedout in a stepwise manner 
assuming a series of differential sections of the evaporator tube and 
using the criteria of a maximum in liquor temperature profile and a 
minima in the pressure gradient curve for the transition from nonboiling 
to nucleate boiling and from nucleate boiling to annular film region 
respectively. 

Simulation programme developed in this si?lcly is compared with 

(a) (2 45) 

the experimental results of Brooks and Badger and Gudmundson ’ for 
evaporating water and subsequently used for the concentration of 
bamboo and baggnse black liquors- 



3 


i -gM.a g. ■ 

' fltiid llfg Battlam f laat iCmnsfeg Meohanlros in a Long 
. Tube Vertioa,! Svapopator 

LTV evaporators generally consist of 2.5 - 7*5 cm. diameter, 

5-10 meter long tubes. Fluid flow patterns and heat transfer mechanisms 
obtainable for a liquor feed below its saturation temperature are shown in 
Pig. 2.1. The following flow patterns and corresponding heat transfer 
meohahiiBs are normally observed* 

1. Single Phase Liquid Flow (Conveotive Heat Transfer); Solution is 
heated to its saturation temperattire by uSonvectlon in the ncnboiling 
section with the cold feed introduced at the bottom of the tube* 

2. Bubble Flow (Nucleate Boiling) ; Superheating of the liquor adjacent 
to the tube wall occurs to some extent and small bubbles of vapor form 

on the heat transfer surface and res’ults in a bubble flow pattern as they 
are released from the surface. 

5. Slug Flow (Nucleate Boiling and Convective Heat Transfer); With an 
increase in number and size, the bubbles coalesce and two phase slug flow 
commences up the tube* In the beginning heat transfer takes place by 
nucleate boiling mechanism but with the increase in slug size, heat transfer 
takes place through liquor film’ by forced convection mechanism*^ 

4. Annular Film Flow (Forced Ccnvection) ; The increase in vapor velocity 
results in a centra,! core of vapor flow and an annular climbing film of 
liquor and heat is transferred across this liquor film by forced convection. 

5, Mist Flow (Forced Convection) * With an increase in temperature gradient, 
dry patches begin to appear on the tube surface and liquid droplets are 
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2,t“ Fluid flow pattern and heat transfer in a 
, LTV evaporator. , ' ■ 
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carried as a mist "by the high velocity vapor. Heat transfer takes place 
through liquor film and at places of dry patches, directly to vapor by 
forced convection. 

6. Single Phase Vapor Flow (Forced Convection); As the vapor rate is 
increased liquor entrainmoat becomes more and more predominant until, 
finally, the liquid annulus disappears and all the liquid is carried in 
the vapor core, giving essentially vapor flow up the tube. 

In a commercial LTV evaporator, a certain length of the tube 
is used in heating the feed to the saturation temperature; in normal 
operation bubble flow, slug flow and annular film flow regions are 
encountered. Mist and single phase vapor flow patterns are not usual in 
chemical process evaporators unless very large temperature gradients are 
impressed on the system; further evaporation in a single pass is limited 
to less than 50~6C9S. Thus it is necessary to consider regions (l) to ( 4 ) 
listed above in a detailed analysis of heat transfer in a LTV evaporator. 

The flow regions and the accompanying heat transfer mechanisms, 
shown in Figure 2.1 are not very distinct and there is considemble overlap 
in actual operations For the purpose of this analysis the following heat 
transfer mechanisms are considered: I) Convection (nonboiling section) 

2 ) Nucleate boiling (bubble flow and part of slug flow region) 5) Forced 
convection (remainder of slug flow and annular film regicsa) • 

Differences in fluid flow pattern and heat transfer mechanisms 
would result in a variation of the heat transfer co-efficient along the 
length of the evaporator tube. There is a sli^t maxima in liquor 
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temperature caused ty pressure drop indicating the start of nucleate 

■foiling region after the nonboiling section. Likewise transition from 

sing flow to the annular film flow regime is accompanied by a minima in 

the pressirre gradient curVe along the tube caused by the high void fraction 

of va,por (0.8-0. 9) which in turn leads to high frictional pressure drop. 

(l5') 

This has been observed by Anderson et.al.' in their experimental 
findings. These criteria may be used to predict approximately the sections 
of the tube lengths whepe single phase liquid flow, bubble & slug flow, 
and annular film flow regions occur during evaporation. . Heat transfer 
co-efficients in these sections can be predicted using the appropriate 
available .correlaticns. These are discussed in Chapter 4 . 

As the liquor traverses up the tube, liquor temperature increases 
due to 'sensible heat transfer and pressure decreases as a resvilt of static 
and frictional forces in the nonboiling section; in this section liquor 
concentration remains constant as there is, no evaporation. Changes in 
liquor tanperature, pressure, pressure gradient and heat f lux in a differential 
length of the nonboiling zone can be predicted by continuity,, momentum & 
en'ergy balance equations? and frictional pressure drop .equation for the 
sin^e phase liquid and heat transfer correlation. The transition to 
nucleate boiling occurs when the liquor temperature is equal to the saturation 
tempera tiure coiresponding to the pressure at tha.t .point in the tube, - In 
the nucleate boiling .regien two pins % flow starts with .an incre.ase in the 
voltime of the vapor by evaporation aoeompanied by an increase in the liquor 
concentration, A differential length in this region can be considered for 
predicting the changes in liquor temperature and concentration, pressure, 
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pressure gradient, quality, void fraction and heat flux by applying 
continuity, momentum & energy balance equations alongwith void fraction, 
two phase frictional pressure drop and heat transfer correlations. 

The stepwise procedure can be repeated until a minima appears in the 
pressure gradient along the length indicating the start of the annular 
film flow region. Changes occuring in this section can likewise be 
predicted by applying continuity, momentum & energy balance eqioations 
along with other equations to differential segments of this region in 
a stepwise manner until the evaporator tube length has been traversed. 
These stepwise calculation procedures are discussed in. Ohapter 5. • 
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CmPTER 3 

SHTERGYt MOMEINTIM & CCSjlTIinjITY BlIANCES gSCTft 'HOITBOIT>ttj<;. 

BOILING REGIONS 

Baergy, Mcanentum & Continuity balance eqtiations are necessary 
for the simulation of the LTV evapcccator. These equations will differ for 
the nonboiling section (sin^e phase) and boiling section (two phase flow) 
of the evaporator tube length. 


The appropriate balances applicable to the single phase region are 
discussed first. In the two phase flow regime energy, momentm & continuity 
balance esqwessions are obtained for (1) homogeneous model and (2) slip model. 
The former model assumes that liquor and vapor travel at same velocity 
where as the latter model is developed on the basis of different velocities 
for liquor and vapor phases. 

A . SINGLE PHI SB LIQ.UID 


The following assumptions are made for the one-dimensional 
analysis of the single-phase flow ewisting in the nonboiling section of a 
e 


A 

dZ 

VK 


(T 

+dT) 

(P+dP) 

! - 

g/gp 

si " 

t 

T 
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27Fr.dZ 
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evaporator ttibei ' 

1 f There is no radial variation of pressure, temperature, Velocity 
& density. ' / 

2 . Change in the Kinetic Energy is neglegihle. 

filer gy Balance? . 

Consider a differential section of the nonboiling section (dZ) 
shown above* Bqnatiajl. (5*1) gives the general differential form of the 
energy equation. 

d(Enthalpy) + d(Kinetic Energy) + d(Potential Energy)>=dq-dW ( 5 * 1 ) 

The terms in equation 3.1 can be evaluated separately as described 

below. 

dH= Enthalpy' of mixture leaving differential section 
- Enthalpy of mixture enetering differential section 
= C,. dT , B'tu/lb-m " ( 5 . 2 ) 

Differential change in Kinetic Energy = 0 ( 3 . 3 ) 

Differential change in Potential Energy = dZ g ( 3 * 4 ) 

J «c 

The amount of heat transfered per pound Hnass of fluid flowing 

is given by equation 3.5* 

Q.2 -fjrr . dZ Q * 2 . 7 T r . dZ 

dq - (3.5) 

w . 5600 f . U , c . 5600 

, h L- ' 

The amount of work done by the system on the surroundings 

dW =0 (5.6) 

S ' ^ 

Substitming equations 3*2 to 5.6 in equation 3 . 1 * energy equation 



gives* 


dT 


Q, 2rr . dZ dZ g 


(5.7) 


Ectmtion 5.7 caji be used to predict the tenperatiore change in 
a differential section of the nonboiling zone. 

Momentum Baianrio t 

Application of Newton's second law of motion to the fluid in 
the differential section (dz) gives the maaentum balance equation 5.8. 

Fotoffl of accelsration - forces of deceleration 

« change in momentum of fluid (5.8) 


as 


Change in the momentum of fluid is 

^ (u + du) - j^u U 


«c 


(5.9) 


When the equation 3,9 is stated mathematically, it becomes 

C^P - (P + dP) c. - ^ ^ 

^ ^ - ^ ^ ^ - 0 (5.10) 


Frictional pressure drop and shear stress qt wall can be relo-ted 




'dp. 




w 


I) 


(3.11) 


P dZ 4 

Py using equation 3*11 , . equation 3. 10 can be rewritten as 


dP 


-g ^ ^ •, 


dZ 


C3.12) 


Equation 3.12 can predict the change in pressure along the tube 
length for the honboiling seotioii#^ 
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For making a one-dimensional analysis of tbe two-phase flow 
in the boiling section of an evaporator tube, the following assumptions 
axe made: 


1. No radial variation of pressure, temperature, velocity, quality and 
specific heat (This implies homogeneous vapor liquid mixture) . 

2. fielqtive velocity between the vapor and liquid ie equal to zero. 

(no slip occurs between the vapor and liquid). 

5. Vapor and liquid are in equilibrium. 

4. Density of vapor-liquid mixture is given as 

f - ( 1 - 0( 



Consider a differential section of the nonboiling section (dZ) 
shown above. The general energy equation 5*1 is applied th this differential 
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section ly evaltiating each term separately. 

dH » Enthalpy of mixture leaving differential section 

Enthalpy of mixture entering differential section (5*15) 

The mixture leaving the section with a tanperattire (T + dl) and a q,uality 
(X + dX) had an enthalpy of . 

( l-X-dX) (T-T^+dT) + (x+dx) ^ + (X+dX) (T+dT-JT^) 

■= (1-X) (T-T„) + (l-X) dT - (T-T J C- dX + X > 

+ dX + X Cp (T-^q) + (T-T^) Cp dX (5.14) 

The mixture entering the section with a temperature T and a qiiality X had 


an enthalpy of 

(1-X) (t-^^) + X^+ X Gp (T-T^) (3.15) 

dH can he obtained from equations 3*14 and 5.15» and ignoring second order 
differentials gives equation 5.l6. 


dH = (1-x) Cp + > dX + X Gp dT 


(if T^-T) 


The differential change in Kinetic energy (dK.E.) is 


d(K.E) 


(u ... 

2 




2 g J 


u du 

®c J 


The differential change in potential energy (dP.F.)i8 


d(P.E.) 


J e. 


(5.16) 


(3.17) 


(5.18) 


The amount of heat transfered per pound mass of fluid flowing (dq) is 


aq 


Q . 2 TTr . dZ 


3600 


W, 


Q . 2 r , dZ 
f u G^ , 5600 


(3.19) 
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TIxg amount of work (dW^) done by the systaoa on the surrounding is 

(3.20) 


dW = 0 
s 


Combining etjuations 3*16 to 3.20 and eq.uation 5*1 gives the energy etjuation 
for twoiphase flow. 

( 1 -x) Cdi + ^vdx + xcp m+ az « . a. 2Trr. az . 

fl d, , 5600 

Elq.mtion 3*21 can be used to predict temperature change or heat flux 
in the differential section. 

Momentum Balance : 

Momentum balance expression represented by eqmtion 3*8 can be 
applied to the differential section to give equation 3.22. 


P^( P + dJP) - f M...g - Y'lr TT D dz 


jf u (u+du) 


fti C, 


u 


S, 


( 3 . 22 ) 


■'C ^ 

Ey using equation 5.11» eqtiation 3*22 can be rewritten for pressiare drop 


in the differential section 


- dP 


f 


■gg„g.. 


g. 


/ ^ ) 

■hf-} 


dZ + 


f u 


du 




(5.23) 


-C ' f '-‘0 

Equation 3 *23 can be used to predict pressure drop in the differential 
secticsi. 

Continuity Balance ; 

Continutiy eqxjation for the boiling section is represented by 
equation 3.24. 

\ f ^ (3 .24) 
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This can be given in differential foim for steady flow conditions in the 
boiling section ty equation 5,25 

n d|^+ j^du =0 (3.25) 

Liquor Solid Concentration ; 


(SfdS) 

(Wj^+dWj.) 


(W^+dWg) 



Solids balance for the differential section is given by equation 3.26 


Total solid entering = Total solid leaving 


S W = S ¥_ + S dW_ + dS + Sy, dS 
^ L L L 


( 3 . 26 ) 

(3.27) 


Change in liquor concentration over the differential length of the 

i 

evaporator is given by equation 3,28,. 


dS 


“S d¥. 


¥T+d¥_ 
li il 


Mass balance for the evapoi?ation gives 


dSL 4- (IaL = 0 

Jb VJT 


Vapor q-uality of the mixture leaving the section is 


X 4* dx *=* 


\ 


a X + 


d.¥. 




(3.28) 


(3.29) 


(5.30) 



Eq.mtion 3 *29 and 3*30 can be combined to give 

+ \) « (5.51) 

By using equation 3»3'' and 3*28, equation 3*32 can be obtained to predict 
the change in liquor concentration over a differential se^ent of the 
boiling section of the tube. 

S(W_ + ¥„) dX 

ds . — - — 2 — (3.32) 

?L- (%+^n) 3X 

iJ IJ (i 

C. TWO PHASE FLW (SLIP MODEL) 

This model differs from the hcanogeneous model in that there 
is slip between liquid and vapor which is taken into account by the slip 
model. The following assumptions are made in the analysis. 

1. No radial variation of pressure & temperature. 

2. Vapor and liquid are in equilibritfei. 

5 . Density of vapor-liquid mixture is given by 

f . ( 1 

Continuity Balance: 

Equation 3*55 cau be written for continuity balance 
W^ = + ¥^ (5 •53) 

If the areas of ccoss section of the two streams are A_ and Ap* 
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The mass flux of each stream is then 

^ W (1~0^) (3.56) 

Gg - P(. vg ot (3.57) 

using the definition of vapor quality (x), the above two equations can 
be used to give alternative expressions for the overall mass flxix. 

(3.38) 




ft? 


X 


Momentum Balance t 

Many alternative forms Cf the momentum equation can be derived 
by manipulating relationships among void fraction, vapor quality, total 
mass fltix, liquor velocity, vapor velocity and other variables. For 

( 53 ) 

steady flow in a pipe Wallis 


dP 

dZ 


-dP^ 


dZ 


”i7 


has given the following expression. 

J + £ 0 ^ f + (1 - ^ ) /j. ] ( 5 * 39 ) 


Combining equation 3.59 with equation 3*38, momentum equation for the 
differential section can be written 


dP, 


dP 

dZ 


P 


G" 


1 


dZ 


dZ 


r 


r' d-x)^ + ^ 

I 7r< j 


(3.40) 


(1 -4) f L + 0^ /g J 

Equation 5f40 can be used to estimate pressure drop along the tube. 


Energy Balance : 

In energy balance, all terms will be alike as in the homogeneous 
model equation 3*21 except kinetic energr term 
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d (K.E) = -i 


X V. 




(1-X) 




(3.41) 


ComlDiniiig equation 3.38 and 5*41 » Kinetic energy term becoraes. 

^ ^ f 0^ 7 


d(K.E) 


X' 


(i-x)- 


•yrn m ^ 


2^0 J 


2 — :?■ 


fS 


(3.42) 


G ‘ fj^ (1-^ J 

Subitituting all terms in equation 5*1 » energy equation becomes 


(l~X) C_dT + XdX + XC dT + S! 

Jj ^ t) rk 




~ LTT^ /l"" 


(1-x)- 
~T 


+ iz 

J e„ 


(k 2 Jf T . dZ 

u . 5600 


(3-45) 


Equation 3.^3 can be used for predicting temperature change 
in the differential section during two phase flow assuming slip model. 

She change in the liquor ccncentration for this model will be the same as 
for the hcmogeneous model and is given by equation 3 ' 32 . 

Equation 3*1 to 3*12 are used for the nonboiling section 
and for the boiling section either equations3.13 to 3*32 assuming homogen- 
ous model for two phase flow or equations 3*32 to 3*43 bas4d on slip model 
for two phase flow can be used. 
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REVIEW OF HE/iT TRANSFER. PRESSIIR.E DROP /\¥D HOLP-HT COHREIATICMS 
'iKD LIQUOR PROPERTIES 


Fluid flow patterns and heat transfer mechanisms in a LTY 
evaporator were discussed earlier in Chapter 2« Changes occur in liquor 
conditions such as temperature, pressure, pressure gradient, concentration 
and factors relating to two— phase flow hehavior. such as quality and void 
fraction, and associated heat flux and overall heat transfer co— efficient 
along the length of the evaporator tube. All these quantities are 
determined in simulating the LTV evaporator. These can be predicted for 
a given evaporator unit from a combination of the following expressions: 
continuity, momentum & energy balance eq\aations, and frictional pressure 
drop and holdup correlations for single phase and two phase flows, heat 
transfer correlations for the nonboiling, nucleate boiling and annular 
film regions, and correlations for the engineering properties of the 
process liquor such as thermal conductivity, viscosity, density, surface 
tension, boiling point rise and specific heat. 

Momentum, energy and continuity relations for one dimensional 
single phase and two phase flow were discussed in detail in Chapter 
This chapter deals with the following aspect l) Critical review of heat 
transfer correlations for an LTV evaporator and recommended expressions 
for heat transfer, 2) frictional pressure drop and hold-up correlations 
to be used in this analysis and 5) Engineering properties for the process 
liquors. 

A. Review of Heat Tri j- nsfer in a LTV Evaporator; 
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include convective heat transfer, nucleats boiling and fohced 
convective heat transfer. A summary of the available literature 
pertaining to these heat transfer mechanisms is given in Table 4.1* 

Available heat transfer correlations are summarised 
separately for each of the three sections - non boiling, nucleate 
boiling and annular film regions. The complete correlating equations 
are not given; instead only the major variables that are included in 
the various investigations of heat transfer in LTV units are given. 

The tube size, liquor systems, accuracy of predictions and the range 
of application of these correlations are shown in Table 4*1 fo^^ the 
various expressions proposed by different investigators. It is 
necessary to choose expressions appropriate for the simialation of 
this investigation from amongst the various equations in Table 4»1» 

The use of empirical overall coefficients representing the cumulative 
effects of heat transfer in the different zones and the steam side 
characteristics are not suitable for the simulation purposes (Equations 
in section A, Table 4»1)» The liquor flow in the non boiling section 
of a natural circulation LTV unit is in the laminar flow region for 
which the non boiling heat transfer correlations proposed by Badger 
et.al.^^^ and Martinelli & Belter can be used. Both of these 
expressions give essentially the same value of the non boiling heat 
transfer coefficient even though the latter include the effect of 
natural convection in their correlation (Appendix B) . For computational 
purposes it is simpler to use Badger et.al. correlation than 
Martinelli & Bolter expression; in the latter expression non boiling 
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section length is an implicit variable which will req.uire a trial & 
error procedure fca? solution. 

The heat transfer correlations given (Section C, Table 4»'l) 
represent the boiling heat tra^nsfer co-efficient for the ccanbined 
nucleate boiling and annular film regions and hence cannot be used 
to predict separately the nucleate boiling & annular film region heat 
transfer rates. 

However , correlations have been proposed for the nucleate boiling 

region by several, investigators (section D, Table 4*1)* Guerrieri & 

Talty^^®^ correlation requires laminar film thidcness which is often 

difficult to determine from opecating units. Correlations proposed 

( 11 ) 

by Dengler & Addams applies to forced circulation vaporization only. 
Mcnelly & Coulson^^^^ expressions cover a wide range of equipment 
geometry & liquor systems and includes the effect of all major para- 
meters. These authors have given two correlations for the nucleate 
boiling region applicable to low heat flux and highest flux conditions, 
and can be usedy';^edict heat transfer coefficient. 

Heat transfer correlations proposed by various authors for 
the annular film region are summarised in (Section E, Tabl2 4*1)« 

The prediction based on Guerrieri & Talty^'^®^ correlation is rather 

f 1 1 ) 

poor and Bengler & Addoms^ ' correlation is applicable to forced 

circulation units only. Anderson et.al.^^"^^ theoretical analysis can 

(19) 

be used but requires a rigorous conputational procedxire. Kuna 
theoretical analysis has not been substantiated by adequate experima^tal 
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fi6') 

data. Penman & Tait^ correlation is based on the use of’tTro 
dimensionless groups and has been shown to fit very well with their 
experimental data and compare favourably with the predictions from 

(14) 

Anderson et.al. theoretical. approach. Hence correlation proposed 
by Penman & Tait can be used to simulate the heat transfer in the 
annular film region of a LTV evaporator tube. 

A summa,ry of the heat transfer correlations based on the 
above review for simulating heat transfer phenomena in a LTV 
evaporator tube is given belows 


1 . Convective Heat Transfer to Single Phase Liquid (Non boiling Hegion^ 

( 6 ) 


Badger et.al, 




/ 0 2. \ 


\ NB 


(4.1) 


U„„ = Overall heat transfer coefficient for non boiling 

NB 

section, BTU/hr st ft. 

C = Specific heat of liquid, BTU/lb “P 

^ Tj^g = Mean temperature difference for non boiling section, °P 
= Density of liquid, Ib-m/cu.ft. 

JJ = Viscosity of the liquid Ib/ft. hr. 

( 10 ) 

2. Nucleate ^oiling Heat Transfer Coefficint- Coulson & Mcnelly 


Low heat flux equations 

. (1.3 . 39d) /_A'f * 


V 


^ 1 i 


Lm 


( 4 - 2 ) 


L / 



27 


N-u = 




i,/i ,fr 


, Re^ = 




Tf ^ 


High heat flux equations 

, 0.69 < 

h^d / / Q. d I 

-— = 0.225 / r j I 1 


P. d ■ 


(4.5) ■ 


^ = 


= Internal dia.,. ft. 

= Nuclaate Boiling heat transfer co-efficient, BTU/hr . sq.ft. °P 
= ThcOTial conductivity i BTIl/hr ft. “P 
= Specific heat, BTU/lh °P 
= Viscosity, Ib/ft. hr 
= Feed rate, Ib/hr 
= Heat flux, BTU/hr sq-ft. 

= Length of the tube, ft. 

= Latent heat of vaporization, BTU/lb 
= Rressure, Ib/sa Tt. 

= Density, Ib/cv ft. 

; } : : — V ' 

= Surface tension, Ib/sq hr 
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5* Annular film heat transfer coefficient - Peamnan & Tait 


^PL 


0.012 


1 0.5 


J 


(4.4) 


Where 

h = film heat transfer co-efficient BTtT/Sec. sq.ft. °P 

= specific heat of liquor, BTU/lb °7 
ph 

D = Diameter, ft. 

J' = Density, Ib-m/cu ft. 

(r~ = Surface tensicsn, Ib-f/ft. 

7^ = Velocity of -vapoiir, ft/sec. 

g = Constant, 52*2 Ib-m. ft/lbf. sec 
c 

4. Shell” Side Heat fransfec Co-efficient- 

The shell side heat transfer co-efficient for the filmvdse 

condensation of steam can be determined from equation 4.5 (Mcadams^^^^) 


h = 1.13 


L Me (Tst- ^3^ 


(4.5) 


Where 


^ H = Latent heat of condensatien, BTU/lb 

2 

g = gravitational acceleration, ft/hr 

L = Length of tube, ft 

Tg '= Surface temparati»ce^°p , . 

= Saturated vapor temperature, "P 
S V 

D = Outside diameter of tube, ft. 
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B. Frictional Pressure Drop & Hold-tip Correlations s 

Correlations representing frictional pressure drop and hold 
up are necessary for the prediction of pressure along the length of the 
evaporator tube length "based upon energy, momentum and continuity balances. 



= 4 log 



2.28 - 4 log 




(4.9) 


For the case of nonisothermal flow in pipes, there is small 

(29) 

change in the value of friction factor as given by Perry' and 
Seider et.al.^^*^^. In the case of heating 


Divide f "by ( for laminar region (4.10) 

a ' w 

Divide f by { M J ^0.17 turbulQit region (4*11) 
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(The units to use in equations 4.6 to 4 .II are as follows: 
D = Tube Dia., ft 

L = Tube length, ft 

= Fluid density, Ib/cu ft 
V = Fluid velocity, ft/sec. 

f = Friction factor 

= Pressure drop due to friction Ib^/sq ft. 
d = Tube dia. inch 


w 


Radius of tube, inch 
Viscosity at mean temp., op 
Viscosity at wall temp., op ) 


The prediction of pressure drop in the boiling section would 
involve a detailed consideration of two-phase liquid-vapor flow patterns. 
Frictional pressure drop & hold-up correlation for the two phase flow have 
been reviewed by several authors^^"'’^^*^^*^^*^^"^”'^ and have reccmmended 
the use of Lochart & Martinelli^^^^ correlation for two phase frictional 
pressure drop given by equation 4*12. 


/ A ^ 


J ntr 1/- / .. V- 


fric 


P 


5-n 


n 


f'lS 


D, 


(4.12) 


The calculation of two phase flow pressure drop by Lochart & 
Martinelli analysis is discussed in Appendix A. 

Levy^^^^ correlation can be used to prectict the hold-up for 
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the flow of air water & liquid-vapor mixtures in vertical pipes. 


( 1-X )' 




( 1-X)' 


1- o( ^ f-Q . 2 ( 1_ 

Where X = Quality of vapor 

= TZhpor voluii).e fraction 

= Density of liquid 

„ = 'Density of vapor 


( 4 . 15 ) 


C. Engineering Properties of Process Liq^or3^ 

The use of momentum, continuity and energy, frictional pressure 
drop, hold-up & heat transfer correlation for the simul.ntion of long tube 
vertical evaporator requires a ccfeipilation of the engineering properties and 
charecteristics of commercial process liquors,. These properties include e 
density, viscosity, specific heat, thermal, conductivity, surface tension, 
boiling point elevation and vapor pressure, data. This study deals with the 
evaporation of black liquors from alkaline pulping processes. Paper industries 
in India use bamboo as the main cellulosic raw material . Other fibcrous 
raw materials used in current practice, include mixed hard--wo6ds, bag^se, 
grasses, straws, salai and pinewood; the use of these raw materials is rather 
limited and is determined by their availability to the paper mills. It is 
hoped that in the next four-five years the use of eucalyptus will be gradually 
established as a important pulping raw material. Several mills use bamboo as 
the mai or ■ raw material and meet the balance of their requirements from same 
or the other cellulosic raw naterials listed above. These raw materials are 
...n,prHially pulped separately owing to diverse differences in their physical. 
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chemical and morphological properties* The spent liquors from these pulping 
operations are usually mixed and treated for chanical recovery. The weak 
black liquor contains 12~2^ dissolved solids approximately two-thirds of 
the black liquor solids are organic compounds dissolved from the 
raw mata:ial during alkaline pulping reactions. The remaining one third 
portion of the black liquor solids constitute the spent inorganic chaaicals. 
The organic compounds include lignin and ca.rbohydrate degradation products, 
and are dissolved in black liquor as their sodium salts. Other orginic 
ccmpounds include soaps (sodium salt of resin and fatty acids present in the 
wood). This often leads to foam problems during handling of -weak black 
liquors. The inorganic compounds in the black liquor include residual 
sodium hydroxide, sodium sulfide, sodium carbonate and smaller amounts of 
sodim sulfate, sodium thiosulfate, sodium sulfite etc. Black liquors 
from the pulping of bamboo, baggase and agricultuiral residues usually 
contain a appreciable concentration of silica as sodium silicate in solution. 
This usually leads to scaling problems in chemical recovery equipments 
requiring frequent shutdowns for maintenance, when the weak black liquor is 
concentrated to 45--50?^ solid concentration in multiple LTV evaporators 
consisting of 5-7 effects. The concentrated black liquor is subsequently 
burnt in the recovery boilers for reclaiming pulping chanicals and to utilize 
calorific value of the orginic compounds in the black liquor solids. 

Even though bamboo is a major raw material in use today, it is 
possible that substantial quantities of bag^se from sugar mill and 
eticalyptus from plantation may become available to the paper industiy. 

Such a situation would then require that a given pulp mill be in a position 
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to handle two or more .major fibrous, raw materials. The chemical recovery 
units must be capable of handling the black: liquors from different fibrpu.?- 
sources. Koorse^ '' has demonstrated the diverge differences in some of the 
engineering properties specially viscosity of bamboo, baggase and ®'‘^C2-l7Ptus 
black liquors. 

¥e can forsee a conditionin the next few years when a particuLar 
paper mill will be. handling bamboo for five to six months and baggase for 
the rest of the year. It would be possible to predict the changes in the 
evaporation plant on switching from bamboo to baggase using the simulation 
programme of this investigation* This would however require ntmerous 
engineering data for bamboo and baggase black liquors. The latter information 
has been experimentally obtained using commerical bamboo and bag^se black 
liquors from Indian paper mills by Koorse^^^^ . Correlations for the engineering 
properties of the liquors are summarized below: 

(38) 

1. Properties of Bamboo Black Liquor - Koorse 

Yiscosityj yU = Exp ( A + B»S + C.S^ + B.S^ ) , cp (4*14) 

Where A = 4*4014 - 0.2248 T +0.0018 T^ 

B = -0.479 + 0.0296 T - 0.00025 T^ 

C = 0.02075 - 0.0015 T + 0.00001 T^ 

D = -0I0OOI75 + 0.00015 T - 0.0000001 T^ 

Density: “ 1.0182 + 0.00779 S - (0.000219 + 0.000003 S) . T, g/c 

(4.15) 

Therml Cc(nductivity:K_ 0. 55354269 - 0.005425 S + (0.0011555 +0.00000209 S)T 

KCal/hr. m^K (4.I6) 
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Boiling point Elevation: ■ B.P.E. = 0*19 3 » ®C (4*1T) 

Eq.uation 4 •''7 for boiling point rise ms obtained from limited 
experimental data and hence needs further refinement with additional data . 
Hevertheless it can be used as an approximate expression to include the effect 
of boiling point rise in the evaporator calculations. 

(58) 

2. Properties of Baggase Black Liquor - Koorse''^ 

2 

Yiscosity; = Exp (A + B.S + G.S ) , cp (4*18) 

where A = 18.276186 - 0. 507619 T + 0.00345 ' 

B = -1 .427425 + 0.045054 T - 0.0003074 

2 

C = 0.052977 - 0.00898 T +• 0.0000064 T 

Density; = 1.0207+0.0071 S-T (0.000279 + 0.0000012 S) g/c6 (4.19) 

Thermal K = 0.050579-0.002563 S +(O.0OO899+0. 00000749) S 
Conductivityt KCal/hr.m°K .. (4.20) 

Specific Heat: 0=0 .995-0.01 258 S +(0.0000556+0.000104 S) (4.2l) 
h 

Boilding Point Elevation: B.P.E. = O.14 S » °C (4.22) 

Equation 4.22 can be used as an approximate expression to include 
the effect of boiling point rise in the evaporator calculations. 

(55) 

3, Properties of Pinewood Black Liquor - Gudmundson 

Viscosity: = Exp ( A+B.S + G.S + D.S ) 1 cp (4*25) 

mere A = 0.4717 - 0.02472 T + 0.7059 x 10"^ T^ 

B ='0.06973 - 0.5452 X 10"^ T + 0,1656 X10~^ T^ 

C = 0.002046 + 0.3185 X 10’'^ T - 0.9761 X ^ 0 ~'^ 

D = 0.5795 XI O"*^ - 0.6129 XIO"^ T + 0.1837 X 10“® T^ 
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Density; f “• 0»495 T + 6.0 S , Kg/m^ (4»24) 

Specific Heat: C_ = (1 .0-0.0054 S) 0.418? , WS/Kg °K (4.25) 

L 

4* Properties of Water - 

Viscosity; = 1.996531 - .0160284 T + .00006745 , cp (4.26) 

-6 

Thermal ConductiTity ; K = 587(1+0.00281 (T-20)) 2.389 X 10 . 

Cal/cm.sec. ®C (4*27) 

5* Vapor Pressxire-Temperatijre Relationship - 

Vapor pressure and temperature can he correlated hy Antoine 
eqmtion 4.26 or 4*27 


log^Q P = A - 

B 

(4.28) 

T-B.P.E. +C 

rp = B 

- C + B.P.E. 

(4.29) 


A - ' 

Values of Antoine equation constants are given in Table C-1 & 
C-2 of Appendix C respectively for equation 4*26 and 4.27* 

(The units to use in equations 4*14 to 4*27 are as follows 
P = pressiire, mm of 'Hg. 

T = Temperat-ure, °C 
S = Solid concentration, vrf'fo) 
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CffiPTER 5 

CCMPmTIOML SCHEME FOR LONG TUBE VERTia^L EYAPORiTOR 

In simtilating the performance of a LTV evaporator, process 
feed liquor conditions such as liquor flow rate, concentration of solids, 
temperature and pressures are usmlly available. Simulation consists in 
predicting the profiles of liquor temperature, concentration, heat flux and 
overall heat transfer co-efficient along the length of the evaporator tube. 
These computations would also give variation in pressure, pressure gradient, 
void fraction and vapor quality along the tube length. Computations would 
begin with the feed inlet conditions at the bottom of the evaporator tube 
and proceed through the nonboiling, nucleate boiling and annular film 
regions, and thus traverse the entire tube length. 

Simulation computations are carried out in-, a stepwise manner 
assuming a series of differential sections of the evapoirator tube and 
using the criteria of a maxima in liquor temperature profile and a minima 
in the pressure gradient curve for the transition from nonboiling to 
nucleate boiling and from nucleate boiling to annular film region respectively* 

Eor a weah black liquor feed below its saturation tempera tirre, 
flow chart of computation for the nonhoiling section is given in Figure 5.1 • 
Flow rate, concentration, temperature and pressure of feed liquor, steam 
chest temperature and tube geometry are fixed for starting the calculations. 

A diff erential length dZ of tube (.15 ni®i'er) is considcared to determine the 
heat transfer rates. The liquor temperature in this differential section 
is assmed to be isotheimial at first • The change in pressure in this section 
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Fig. Flow Dia^cam fear Wonbollirxs Section of a JEY Tube 


i START 

i _ ± “3Z1_ 

m,,TS,TI,S,P 


V R^d, XLjDZjTP 


J_ Take a differential step of length dZ 
T 3DT = 0 


Compute frictional pressure drop 

Eq. ( 4 * 6 - 4 . 11) 



Compute total pressure drop by 

momentum Eq.(5.12) 


Compute overall heat transfer 
co-efficient. Ea.(4*l) 


I Compute Heat~ f lux by Eg,. Q=I JA/nT | 

rpompute new TT by energy Eg* {5>T) 

. — 

Z' Toes TT differ from estimate by ^ 

V more than fillowable (.1°F) z 



Use new TT 


/ Bfint output! variables for the 
differential section- 


"fes feed reached its satvirati 
temperature 


ion \ 
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aan be calculated frcm the frictionaJL pressure drop correlation equations 
( 4.6 - 4*1l) and momentum balance equation. ( 3 -12) . The overall heat 
transfer co-efficient for nonboiling section can be obtained by iEfedger 
correlation equation 4*1 • He;?.t flux can be computed by the hait transfer 
rate eq'uation ( Q = UAa.T). The next step consists of determining 
differential change in the temperature by energy equation 5 *7* This new 
Value of DT is then compared with the previously assumed value for the 
start of the computations. If they differe by more than a certain 
specified limit (.1 ®F), then the calculations are repeated with this nexv 
Value of DT. In case the difference is within specified limit then the 
calculo,ted liquor temperature is compared with the saturation tanpenture 
corresponding to the estimated pressure at the exit of the differential 
secti'On. If the liquor tempendure Ixippens to be below the saturation 
temperature, the ailculations are repeated for the next differential section. 
This iterative procedure is repeated for several differential sections 
until the liquor tonperature equads the saturation temperature of liquor. 
This point corresponds to the end of the nonboiling regime and subsequent 
calculations are carried out for the boiling section. 

For the boiling section, sim.ulation ccanputations can be 
performed by assuming either homogeneous model or slip model for the trro 
phase Vapor-liquid mixture flow as given below: 

1. H^xnogeneous Model : 

Homogeneous model assume® that the velocil^y of liquid and 
vapor are same and there is no slip between liquid and vaper. The 
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computationad scheae for estimating changes in the process liquor conditions 

and heat transfer rate is shown, in Fig. 5 *2 for the boiling section. Conditions 

of the liquor feed to the boiling section would be the mine as exit liquor 

conditions from the nonboiling section. Now consider a differential length , 

dZ ( 0.15 meter) in the evaporator tube length. The calculation for this 

differential section must give the change in liquor temperature, pressure, 

pressure gradient, concentration, quality, void fraction, heat flux and 

overall heat transfer co-efficients. Firstly values for the change in the 

quality, temperature, solid concentration and velocity ( DX,DT,DS,DV 

respectively) for this differential segaent are assumed. Two phase flow- 

frictional pressure drop is calculated by Lockhart and Maxtinelli correlation 

(equation 4.12).' The void fraction at the inlet and exit of the differential 

zone can be cadculated by Levy correlation (equation 4*15) • Obviously the 

be 

void fraction at 'the entrance to the boiling section will/zero and gradually 
increases as boiling progresses along the evaporator tube. The total pressure 
drop' can be obtained by momentum balance(eq.3 *23) • The temperature 
corresponding to the calculated pressure can be determined from the Antoinl 
eq. 4*17 and this would give new value of change in temperature for the 
differential section. Heat flux for the differential section can be 
obtained by the energy equation 3 • 2 '^ * Overall heat transfer co-efficients 
can be calculated ly using eq.4.2,4.5 and 4*5 • New value of heat flux can 
be calculated by’ heat transfer rate equation (q = UA/^T) . This value of 
heat flux can be used in energjr eq. 3*21 to calculate new value of change 
in quality for the differential section and is compared -with the pre-viously 
assumed value of DX. If the difference between these two values of LX is 
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Flow Diagram for Boiling Section (Homogeneous Model ’Pf-P- LTV Tube 
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larger than a certain specified limit (.0005) then an average of these 
two valnesof quality is used and the calculations are repeated from the 
step computing he^^t fliox from energy eq. 5*21 until the difference between 
calculated and the assumed values of DX is within specified limit. In case 
difference ishvithin specified limit, solid concentration change is obtained 
by eq,. 3*52. How value of change in the velocity of homogeneous phase can 
be calculated by continuity eq. 3*25 and can be compared \=dth the previously 
assTmed values of DV. If the difference between these two values is larger 
than a certain specified limit(.0005 meter/sec) tlan calculations are repeated 
with this new Value. Incase difference is lathin specified limit, pressure 
gradient curve is checked for the minima. Index number which will make change 
from nucleate boiling to annular film flow vdll be changed with the occurence 
of minima in the pressure gradient curve. This iterative procedure is carried 
out for several differential sections until whole of the length is traversed. 

For annular flow region, computational procedure is same as for 
nucleate boiling region except Penman and Tait eq. 4*4 of annular film heat 
transfer CO- efficient is used in place of eq. 4.2 & 4.3. 

Slip Model ; 

Slip model implies a different velocity for liquor and vapor phases. 
Plow chart for the simulation of heat transfer and other phenomena based on 
slip model is given in Pig. 5. 5. This scheme of cxanputation based On slip 
model is essentially the same as for homogeneous model except with a slight 
difference that eq. 3*40 and 5.45 are used to estimate pressure drop and 
tempcsrature change respectiv^y for the differential section. 
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?lg. 5»3 Flow Diagram for Boiling; Section(51ip Model) of a LTV Tube 
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RESULTS & JISgJSSION 

The ccmputational scheme outlined in chapter 5 is used to 
simulate the long tube vertical evaporator for the following liquor 
systems! water, bamboo black liquor, ISaggaee black liquor and pinewood 
black liquor. 

Single Effect LTV Evaporation Simulation ! 

A . Application of Simulation Models to Badger* s Lata for feteri 

The accuracy of the simulation programme for sin^e effect 
evaporator calcialations has been compared with the experimental data 
of Brooks and Badges^^^ for the evaporation of water in a LTIT unit 
consisting 6.1 meter long and 5»08/4.45 Cm diameter tube. Brooks and 
Badges' s experimoital data for water are summarized in Table 6,1. The 
variation in process variables along the tube length is determined by 
the simulation programme based on flow charts Fig. 5.1 and Fig, 5, 2 for 
ncwiboiling and boiling section (Hcmogeneous model) respectively. The 
laagth of the boiling section and the corresponding mtes of evapcnraticn 
and heat transfer, and the average overgill boiling heat transfer 
co-efficient (At-g U = ^ \ y/ S ^ ^ obtained from the 

simulation programme givoa in feble 6.2 for the experiments reported 
by Brooks and Badger. Brooks and Badger's results are also included in 
Table 6.2 for comparison purposes. The range and average percent deviation 
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Table 6.2 -• Simulation Results for Hfeter Evaporation Using 
Homofieneous Model 


Data Boiling Section 

Length (jneter) 

Evaporation Ifeite 
(Kg/hr) 

Heat Transfer 
rate in the 
boili-ng Sect- 
ion (KW) 

Average overall 

H^t transfer 
Co-efficient in 
Boiling Section 
(KW/m^^K) 


Computed ■value 

4.84 

39.6 

24.20 

3.04 

1 

Bxperime'h-fcal -value 

5.65 

75.0 

47.80 

3.26 


ia Deviatton 

-14.3 

-47.2 

-49.5 

-7.1 


Computed -value 

4.21 

73.2 

20.13 

3.16 

2 

Experimen'fcal -value 

5.19 

43.4 

27.62 

3.45 


io De-vlation, 

- 18,8 

69.5 

-26.6 

-8.7 


Computed -value 

4 i 28 

52.8 

17*62 

3.07 

5 

Experimaatal value 

5.25 

38.8 

24.40 

5.57 


io De-vdatibn 

• 18*6 

37.8 

- 28.6 

- 10.8 


Computed ■mlue 

4.67 

105.8 

26.00 

3*07 

4 

Experimen-tal -val-iid 

5.48 

7.5.4 

• 47.55 

3.49 


io De-viation 

-15*0 

43.8 

-45.0 

- 11*9 


Computed -value 

3*59 

58.6 

16.90 

3.14 

5 

Experimen-bal -value 

'5.45 

29.3 

18.10 

3.15 


io De-viation 

- 19.6 

101.0 

1 

• 

0 

0.5 ; \ ^ ^ 



46 


(T^ble 6.2 continued) 


Computed "value 

4.04 

91.5 

22.94 

3.15 

6 Experimental -value 

4.91 

59.5 

57.50 

3.34 

io Deviation 

-19*5 

54.0 

-59.1 

-5.6 

Ccmputed value 

4.25 

79.5 

21 .44 

3.17 

7 Experimental value 

5.12 

54.6 

34 

5.25 

^ De-viaticn 

-17.1 

46.5 

-38.5 

-2.54 

Computed Value 

4.58 

74.2 

21.70 

3.14 

8 Experimaital value 

5.45 

56.5 

55.60 

3.24 

io Deviation 

-16.2 

52.0 

-39.8 

-2.8 

De"v±ation range i 

14-20 

52-101 

5-50 

0.5-12 

io Average De-viaticn 

17.2 

53-9 

53.9 

6.2 


At- 

of the calculated values from the Brocks and Badger’s expsrimental 


observations are given 

below (assuming homogeneous 

model for the boiling 

section) : 



Process Variables 

Deviation range, io 

.Average Deviation ! 

Boiling section length 

14-20 

17 

Bvapo3ra.tion rate 

52-101 

54 

Heat Transfer rate 

5-50 

54 

Boiling average overall 


heat transfer co-efficient 2-12 

6 


The above results show that while this simulation model can 
predict the overall heat transfer co-efficient quite well, accuracy for 
predicting the boiling length can be considered moderate and the estimates 
of the rate of evaporation and heat transfer by the model are poor. The 
rather large deviations observed may be attributed to the limitations 
imposed by the homogenous model which assumes equal velocity for the 
and liquid phases during two phase flow. 

The slip model for the boiling section was then used to 
simulate the process conditions of Brook & mdger' s experiments based 
on flow chart 5.1 & 5*3- Results as obtained by using this model and Badger's 
experimental reaults axe given in Table 6.3 ,bich also include deviation of 
the predicted values from . the experimental values. The following 
observations can be nade from the results in Table 6.5, concerning the 
accuracy of the simulation programme based on slip model fcr two phase 
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-Unable 6,3 Simijlation results for Tlfe.ter Evaporation. Using 
Slip Model; 


Data 

Boiling Section 
Length (Meter) 

Emporation 

I^te 

(Kg/hr) 

Heat Transfer 
rate in the 
boiling Sect- 
tion(KW) 

Average overall 
Heat transfer 
Go-efficdait in 
Boiling Section 
(KW/m 2 ‘^K) 

Computed mlue 

5.42 

68.3 

43.90 

2.68 

1 Experimental mlue 

5.65 

75.0 

47.80 

3.26 

io Deviation 

-4.1 

-8.5 

- 8.6 

-17.5 

Conputed value 

4.21 

40.0 

25.77 

2.76 

2 Experimental mlue 

5.19 

45.4 

27.62 

3.45 

io Deviati on 

- 18.8 

-7.6 

- 6 .8 

- 20.0 

Computed mlue 

4.28 

35.0 

23.00 

2.78 

3 Experimoatal mlue 

5.25 

38.8 

24.40 

3.57 

io l&viation 

- 18.6 

-9.4 

-6.9 

- 21.0 

Canputed mlue 

4.99 

66.5 

42.40 

2.71 

4 Experimental mlue 

4.58 

75.4 

47.55 

3.49 

io Deviation 

-7.3 

-8.9 

- 11 .1 

- 22.2 

Computed mlue 

5.59 

53.4 

21.20 

2,76 

5 Experimental value 

5.45 

29.3 

18.10 

3.15 

io Deviation 

-19.6 

14.8 

16.4 

- 1 1 .6 
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(Table 6.5 continued) 


Computed value 

4.04 

48.8 

50.94 

2.72 

6 Experimaital value 

4.91 

59.5 

57.50 

5.54 

io Deviation 

-19.5 

-17.5 

- 18.0 

- 18.3 

Computed Yalue 

4.25 

44.1 

28 .40 

2.79 

7 Experimental value 

5.12 

54.6 

54.66 

5.25 

io Deviation 

- 17.1 

-I8.9 

-16.9 

-14.1 

Computed value 

4.58 

46.0 

28.95 

2.78 

8 Experimental value 

5.45 

56.5 

35.60 

5.24 

i Deviation 

-16.2 

- 18 . 4 

- 20.2 

-14.2 

— .y — 

Deviation range i 

4 “ 20 

8 - 19 

6 - 2 ' 

10 - 22 

io Av^erage Deviation 

15 * 2 

15.2 

15.1 

17.2 
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flow: 


Process Variables 

Beviation range, io 

Average Beviation, io 

Boiling section length 

4 - 20 

15.2 

Evaporation rate 

8 - 19 

15.2 

Heat Transfer rate 

6 - 20 

15*1 

Boiling average overall 

heat transfer co-efficient 

10 - 22 

17.2 


The above resxilts show ttat the simulation, programme of this 
study assming slip model can. predict Brooks and ledger’s experimental 
conditions with an overall accuracy cf about 15%* Hence the simulation 
programme based on flow charts 5*1 and 5»5 for the nonboiling and boiling 
section (slip model for two phas e flow) respectively Q&a for the 

subsequent calculations of this investi^tion for emporation of bamboo, 
baggase and pinewood black liquor fran alkaline pulping processes. 

B. Compari s on with Gudmundson's Experimental curves; 

The simulation programme is also used for ccanparing the predicted 
values with the experimental results of Gudmundson^^*^^^^ . The evapomtor 
used in the latter work consisted of four 10 m long stainless steel tubes 
with diameters of 51/46 mm. Gudmundson’ s results for water at a feed rate of 18 
18 Kg/an^ are used in this comparison.. 

Figures 6.1, 6.2 and 6.5 give the predicted temperature profile 
and overall heat transfer Co-efficient along the tubelength for the 
various inlet feed tonperatures of T6°C, 81.5°C and 85°Cj and correspond- 
ing steam temperatures of f'1.5°C, 95°C and 86.5°C. The simtlation 
calculations are represented iy the continuous curves in the above figures 
assuming slip model for the two phase flow region. The predicted tempera ture 










Experimentat 

Predicted 

Steam temp. =91-5°C 
Tube dia.= 51/46mm 
Feed water rate =18 Kg/ 


U(KW/m^ K) 








4?.§y^W' 




Experimenta 

Predicied 

Steam temp.=95“C 
Tube dia.= 51/46mm 
Feed water rate = 


U(KW/m^°K) 


on of predicted values with Gu 
jxperimental data for water. 













Experimental 

Predicted 

Steam temp-sES-B^C 
Tube d)a. = 51/46mm 
Feed water rate s 
18 Kg / Srn^ 


Comparison of predicted values 










54 


profiles agree with, the experimental data quite well; the agreement is 
within 1°C for the two sets of data represented, in figures 6.1 and. 6*3 
and. is within 2°C for the data in Figures 6.2. Experimental curves for 
the over all heat transfer co-efficient diviates somewlat from the 
predicted profiles based on this study. Experimental results in general 
show an average deviation of about 255^ from the predicted results. 

The above comparisons for the heat transfer of LTF evaporator show 
that the simulation model of this investigation can be used to predict 
reliable values of heat transfer rates. Hence the simulation programme 
based on flow charts Fig. 5*1 and 5*5 for the nonboiling and bofling 
section (slip model for two phase flow) respectively are used for the 
subsequent calculations of this investigation for evaporation of bamboo, 


ba'ggaae and pinewood black' liquors. 

C. Evaporation of Pinewood Black Liquor ; 

Weak blade liquor from the Kraft pulping of pine is assimed 

to contain 12 . 05 ^ dissolved solids with a flow rate of 100 tons/hr at a 

temperature of 50°C. The LTV evaporator unit is arbitrarily assumed to 

consist of 300 tubes, 5.08/ 4'.4-5.ciandi^. and 7.6 meters long. Steamc&eet 

pressure is taken to be 306 KH/sq.m. Pressure of the liquor feed is 

assumed to be 104 Of/sq m. A print out of the simulation programme for 

the evaporation of pinewood blado liquor is given in Appendix L-l for the 

res presentation of homogenous modd. describing two phase flow pattern . 

(35) 


Physical property correlations reported by Cudmundson 
in this analysis. 


are used 
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Tbe simulation prcg«»ns gives the output variables ih the 
following order- vapor quality, void fraction, pressure gradamt, pressure, 

temperature, heat flux, overall heat transfer co-eifioient and solid 

• -hnrvnc, Of these variables along the length of the 
concentrebion. The variations of these var 

evaporator tube are shown in Pig. 6.-4 - 6 .<'respectively. It can he 

Observed that of the total length of the evaporator tube (7 . me er , a 

j-v, -pricri Tiauor to its boiling point in the 
2 meter is used in raising the feed liquor 

non-boiling region. 

Kg. 6.4 shows the variation in vapor quality (mass fmctron of 

. j_ .PT cinns the tube length. It 
. 4.1, +rv4-nT vninor -liQuid mixtuce flow) dong m 

vapor in the total vapor nsL d 

can be seen that the mass fraction of vapor increases from ^ 

start Of the boiling section to 0.55 and 0.22 towards the an c 
ass^ing hcmogeneous and slip models respectively, .he difference 
Observed between t. two models in .used hy the oontrihuticn o th^ 

,q,.bio ehergy terms (eq.5.43). - -e hcmogenecus model h.et . 

rkVincip velocity whereas the slxp 
IS based on the use of hc^ogehecnis phase 

tabes into account the contribution to kinetic energy a 

vapor and liquid velocities during two phase flow. 

4-11 r-hnnaes iu void fraction along the tube 

Pig. 6-5 gives the changes ru vux 

Ti +Vnt the start of nucleate boiling section is 
and it can be seen that the 

do Tn void fraction caused by rapid 
aocompained hy a steep in..s » v 

uuol.t.» and bubble ^ .t the end 
essentially the same value for void fraction. 
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of the evaporator tube is in the range of 0.95 ” 0,97* ‘ 

Qhange in the pressure gradient along the tube laagth is 
^ven in Fig. 6,6 . Pressure gradient decreases slightly in the nonboiling 
section and subsequently rapidly in the boiling section. The observed 
decr^se in the pressure gradimt is due to static and frictional pressure 

• i 

drop effects. The lower -value of the pressoire gradient for the homogeneous 
model can be attributed to the fact that in homogeneous model, changes in 
the momentum of fluid is based on the use of homogeneous phase velocity 
whereas the slip model takes into account the -vapor and liquid velocities durin.|: 
during two phase flow to calculate change in the momentum of fluid. In 
the upper portion of the tube pressure gradient becomes substantially 
constant for the homogene-oas model whereas it passes through a minima in 
the slip model which can be attributed to incr^ising frictional pressure j 

drop in the annular film region. This minima in the pressure gradient ' 1 

ccrresponcPto transition from nucleate boiling to annular film region. | 

I 

The pressure gradient curve for the slip model is similar to experimental | 

( 15 ) I 

results reported by Anderson et.al . ; 

Fig. 6,7 shows a gradual decrease in pressure and is caused by | 
static and frictional pressxzre drop contributions. The steady decline in 
pressure in the nonboiling section is due to static pressure drop. In the ■ 

hoiling section pressure decreases at a lower rate due to combined effects i 

of hydrostatic head and frictional losses. The curves in Fig. 6.7 are i 

( 15 ) i 

similar to es^erimen-fcal curves reported by Anderson et.al. . 
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Pig. 6.j8 shows variation of liquor temperature in the tube 
length. It can be seen thit less than one third of the evaporator tube 
length is used in raising the feed frcm 50®C to 102.5°C (saturation 
temperature corresponding to the pressure in the tube) in the nonboiling 
section. In the boiling section eacbending over more than two thirds of the 
evaporator length. The homogeneous model and slip model for the two phase 
flow predict liquor tmperatxn?e within 5°C of each other. 

Hsit flux along the tube length is given in Pig. 6. '9 and the 
values predicted by the homogeneous and slip model are nearly equal, ^he 
decrease in the value of heat flux for the nonboiling section is caused 
by the decrease in temperature driving force as the liquor is being heated 
to its boiling point. The increasing trend in heat flux over the boiling 
section in due to the increase in overall heat transfer co-efficient for 
the nucleate boiling and annular film regions. The transition from nucleate 
boiling to annular film region occures at 4*3 meter and 5 *5 meters from the 
results of heat flux for the homogeneous and slip models assumptions 
respectively. 

Pig.6 .-|0 pdves the variation of overall heat transfer co-efficiant 
along the tube length and shows clearly the three dominant h^t transfer 
mechanisms in a LTV evaporator tube. Out of a total of 7.6 meters long 
tube the nonboiling, nucleate boiling and annular film regions ex';end over 
2.0, 2.3 and 3»3 meters respectively for the homogeneous model and the 
Gcrresponding lengths for the slip model are 2*0,3*5»2.1 meters. The 
transition from nonboiling to nucleq-te boiling and nucleate boiling 



to annular film regions are accompanied by a pronounced increase in the 
Value of overall heat transfer co-efficient. In general the homogeneous 
model predicts higher value for overall heat transfer co-efficient ccmpared 
to the value based on the slip model. It can be attributed to the higher 
value of vapor q-uality for the homogeneous model which results in hj.gher 
vapor velocity and this leads to hi^er heat taiansfer co-efficient than the 
slip model (eq.4.2,4*5& 4*4) * 

Fig. 6.1)1 gives the change in liquor concentration along the 
tube length and predicts an increase from a initial concentration of 
to 18.4^ and 15*4/^ at the tube exit end for' the homogeneous and slip 
models respectively. The higher concentrations observed for the homogeneous 
model is due to the higher value of overall heat transfer co-efficient 
obtained in Fig.6 .‘lOcompared to slip model. 

H • E vaporation of Bamboo Black Liquors 

The simulation programme following the flow charts in Pig. 5*1 
and 5*5 based on the slip model for two phase flow is used for the 
evaporation of bamboo black liquor in a single effect LTV evaporator. 
Physical property correlations reported by Koorse^ are used in these 
calculations. Evaporator dimensions and process conditions are chosen to 
be the same as that used for pinewood black liquor in section C earlier. 

Two caases of evaporation are assumed for the calculations namely (l) 
clean tube (2) tube ■with a dirt resistance of 0.45 “K/EW. The 
dependence of the variable on tube length is given in Fig.6 .12 -6 .1# for 
vapor quality, void fiaction, pressure gcadient, pressure, temperature, 
heat flux overall heat transfer co-efficient and solid concentration. 
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The comparison of the curves 6 . 12^,19 with those in Fig. 6 . 4 - 6 .t 1 

shows tint the natxnre of the curves are very similar for the evaporation of 

Bamhoo and pine hlack liquors. It can he seen from Fig.6 that t*hc 

contrihution of the dirt scale resistance is an important factor and 

significantly reduces the value of solid concentration obtained for the 

tubes with clean heat transfer surfaces. Dirt scale resistance increases 

the tube length corresponding to nonboiling region while decreasing the 

tube length available for boiling region. For the process conditions 

assumed for the evaporation of blade liquor, the nonboiling section 

prevailing over 1 meter length of clean tube increases to 2 meter when dirt 

2 

scale resistance is assumed to O.45 m °K/kw. 

Fig. 6.12 Pvgff ^variation of vapor qualiiy with tube length and 
it can be seen that the q.x^iality increases from zero at the start of the 
boiling section to 0.55 and 0.23 towards the end of the tube for clean and 
fouled surfaces respectively. The'higjier vapor quality obtainable for the 
clean surfaces suggest that the process liquor can be evaporated to a higher 
solid concentration for the given length of the tube. Fouled sxarface lave 
the effect of decreasing the capacily obtainable using ol^ner evaporator 
tubes. 

The graph obtained for the void fraction in Fig. 6 , 1'5 for bamboo 
black liquor is essentially identical to the curves obtained for pine black 
liquor in Fig. 6 .§ assuming equal diirt resistance in the tube. 

Fig. 6.14 gives the variation of pressure gradient along the 
evaporator tube. For the last' part of the tube length value of pressxire 








BAMBOO BLACK LIQUOR 























>Vvr 
















gi®?S?i’.s?jfjf$* 

'C'r'iy-'' 






.<t 

o 

o 


if 

tl 

00 

'•O' 

'TD 

00 

or '' 

or 

CD 

o 

a 

2 


: < 



00 











If 

gradient is higher for the clean tube in comiajsion to fouled tube. The 
curves go through a minima at 2.2 and 5*2 meters for the clean and fouled 
tubes respectively and these lengths represent the transition frcaa nucleate 
boiling to annular film regions. The pressure gradient in the annular film 
region is up to 5 times the value for the fouled tubes having the dirt 
resistance of 0.45 m °K/KW. Such a condition can cause regicns of low 
absolute pressure in the annular film region which may l^d to plugging of 
tubes. 

Pig. 6,15 describes the longitudinal variaticm of pressure along 
the tube length based on pressure gradient given in Pig. 6.14* 

Pig, 6.1$ gives the profile of liquor temperature and shows 
that the maximum value of 106.0 °C for the clean tube and 102.8 °G for 
the fouled tubes. These profiles are similar to these obtained for pine 
black liquor in Pig. 6.$. 

Pig. 6.17 shows thit with clean siirfaces it is possible to 
obtain heat flux rates twice those of the fouled tubes with the assumed 
dirt resistance value. 

The value of the overall heat transfer co-efficient decreases 

2 . 

by a factor of 2 to 5 with the assumed value of dirt resistance 0,45 m °K/K?f 
as demonstrated by the graph in Pig. 6 . 18'. 

The change in bamboo black liquor concentration in Pig. 6.1^ 
shows that a final con centmtiOn of 27 *8^^ and 15»7?S can be obtained for 
the clean and feruled tubes respectively. This points out the need for 
relatively clean tube surfaces so tlat satisfactory evaporator capacity 
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levels ca,n be maintained* 

B« Evaporation of Bag@.se Black Liqnor and a Comparison with toaboo Black 
Liqnor Evaporation: 

The canputations for the evaporation of bag^se black liquors 
in the single effect LTV unit are identical to the procedure used earlier 
for pine and bamboo black liquors. The listing of ccanputational procedure 
for baggase black liquor is given in Appendix I)-2 and is based on the same 
equipment geometry and process feed liquor conditions as used for the other 
black liquors in this investi^tion. Fig. 6. 20-6. 2^ gives the profiles of 
temperature, overall heat transfer co-efficient and solid concentration 
along the tube length for both clean and fouled tube siorfaces. The curves 
marked C-2 and 'D-2 denote clean and fouled evaporator tube (Rjj =0.45 °k/kw) 
respectively for baggase black liquor. For comparison purposes the 
corresponding results for bamboo black liquor are reproduced in these 
figures and the curves marked C-1 and D-1 represent clean and dirty tubes 
(r^= 0.45 m^ °K/K!W) respectively. Baggase and Bamboo black liquors are 
very good examples for danonstrating the role of engineering propo^ties of 
process liquors on the performance of evaporators. 

The curve in Pig.6.20 show that the temperature profile (C-1) 
and (D-1) for bamboo black liquor j and (C-2) and (d- 2) for baggase black 
liquor are similar. The loigth of the nonboiling region is 2.15 meter and 
2.75 meter for the clean and fouled tubes respectively handling baggase 
black liquor and the corresponding value for bamboo black liquor are 0.92 
and 1.85 meter. This shows that the nonboiling section length is increased 
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tiy a factor of 2«5 and 1.5 fo^ clean end fouled surfaces respectively 
when changing over from bamboo to hag^se black liquor. Th5-s can be 
attributed to the lower value of the overall heat transfer co-efficients 
for baggase blade liquor as a result of its higher viscosity. 

The curves in Pig. 6 .21' show that the value of ovarall heat 
transfer co-efficient for bamboo black liquor in the nonboiling section, 
is 2.5 times the value for baggase black liquor. This effect of liquor 
■viscosity on the overall heat transfer co-efficisnt decreases the 
evaporation rate as may be observed by the curviis in Pig.6.^2 (Curves C-1 j (;”2 
and D-1 , D-2) . The viscosity of the Bamboo and baggase black liquor for the 
concentration range of Pig.6.22 are given below in Table 6.4. 

Table 6.4s Pisoosity Bangc for Bamboo & Bxggase Black Liquors in a Clean Tu be 

aimboo Black Liquor Bag^se Blade Liquor 
Inlet Exit Inlet Exit 

1. Honboiling section- 


Tempera turej °C 

50.0 

106.0 

50.0 

102.4 

Solid Concentration j 

12.0 

12.0 

12.0 

12,0 

Tisccsity, op 

0.9 

0.45 

4.7 

1.9 

Boiling section- 





Temperature, °C 

106.0 

95.5 

102.4 

36 .2 

Solid Concentration, io 

12.0 

27.8 

12.0 

25-9 

Viscosity, cp 

G.45 

1 .1 

1.9 

9.8 
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The restilts of the above canparison clearly shows that a 
given evaporator can he effectively used using preheated haggase black liquor 
to reduce tte length of the nonboiling section, ihe role of wide variation 
in physical properties of black liquor can be appreciated in evaluating the 
perfonaance of commercial plants utilizing mtiltiple effect evapcrators 
including auxilary h^t recovery equipments such as liquid preh^ter, heat 
exchangers, condensate vaporizer^ product vaporizers, vapor bleed etc. 
Evaporation plant in the chemical recovery system of kraft pulp mill usually 
consist of 4“7 effects with many of the heat recovery features listed above. 
The ccmputational scheme for the single effect evaporation developed in 
this investigation can be applied to multiple effect evaporators which 
can be suitably modified to incorporate heat recovery features in the 
simulation programme. 
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SUMMARY 

Simiilatian progminme developed in this work based on slip model for 
two phase flow can give reliable values of the heat transfer rates in the 
long tube vertical evaporators. The model has been applied to black liquors 
over the concentration range of 12-28^ dissolved solids and the values of 
the overall heat transfer co-efficients obtained agree very well with the 
experimental values reported in literature. The profiles of liquor temper- 
ature, concentration, heat flux and overall heat transfer co-efficients 
along the evaporator tube length are obtained for bamboo, bagga.se and 
pine black liquors and the influence of solution viscosity on the lengths 

of nonboiling, nucleate boiling and annular film regions in a giveai 
evaporation unit is illustrated. Reduction in evaporator capacity resulting^_ 

from the dirt resistance to h^t transfer from scale deposits can be 

predicted by this programme. The model can also be used to predict changes 

in evaporation capacity and the necessary auxiliary heat enchange equipment 

when the process liquor exhibits a large change in the viscosity. 
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Cross section of liquid stream, Sq ft. v 
Cross section of "vapor stream. Sq.ft. 

Cross sectional area of tute, Sq ft 

Specific heat of liquid, Btu/lb °P 

Specific heat cf vapor at constant pressure, Btu/lh 

Tube diameter, ft 

2 

Gravitatdonal Acceleration, f t/sec 

2 

Conversion factor, 32.2 Ib-ft/lb^. sec 
Mass flux of vapor stream, Ib/sq ft. sec. 

Mass flux of liquid stream, Ib/sq ft. sec. 

Mass flux of liquid-vapor stream, Ib/sq ft. sec 
Mechanical equivalent of heat, 778.17 ft-lb.^/Btu 
Pressure, Ib^/sq ft 

Total pressure drop per unit lengh, Ih^/cn ft. 

Frictional pressure drop per unit length, iP^/cu ft 
Heat flux, Btu/hr. sq ft. 

Solid concentration, Wt ^ 

Temperature, "^F 
Reference tonperature, ®F 
Velocity of hcmogenecus phase 

Velocity of vapor stream, ft /sec 



Vt 




w. 


X 

z 

fL 

fs 

--r 

I w 


A 


Yelocity of liquid stream, ft/sec. 

Mass flowrate of vapor, Ib/sec 
Mass flot? rate of liquid, Ib/sec 
Total mass flow rate, Ib/sec 

Work done by the system on surrounding, Btu/lb-m 
Quality of -vapor-liquor mixture = 

Tube length, ft 

Density of liquid, Ibm/cu ft. 

Denisty of vapor, Ibm/ cu ft. 

Density of vapor-liquid mixture, Ibm/ cu ft. 
Shear stress at the wall, Ib^/sq ft. 

Yoid fraction 

Heat of vaporization, Btu/lbm 
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APPuroix A 

m) PmSE FRICTIONAL PBESSPRE DROP 

(24) 

Lodchart & llartinelli ^ have analyzed four types of two 
phase flow mechanisms which exist during the simultaneous flow of a liq^uid 
and a or vapors (1) 'Turbulent-'Pur'bulent flow (2) Yiscous-turhifLent flow 
(5) ^urhulent-viscous flow (4) ^scous-Viscous flow (flow of both liquid 
and ^s may be viscous) . This analysis is based on the following assumptionss 


1« Static pressure drop for the liquid phase must equal the static pressure 
drop for the ^seous phase regardless of the flow pattern, as long as 
aJi appreciable radial static pressure difference does not exist. 

2. The volume occupied by the liquid plus the volme occupied by the ^s 
at any instance must eqxial the total volme of the pipe. 

3. Flow mechanism transition criteria: 


^ 2000 
gP 

Re, 1000 
Ip ^ 

where gg 

gP 


the ^3 phase is turbulent 
flow of liquid is viscous 



Lock^'"t & Martinelli defined a parameter X (two phase flow 
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vt 


^ -0 
0 , Re 

1 gP 

.8 

1 ' g 

M, 

(A- 2 ) 


C 

g 


¥ r, 
g ^ 1 




C Re 

1 Ip 



f' M 

(A- 5 ) 

^tv 


C 

g 


¥ 

g 

P, A 

1 ^ g 

2 




^1 


(A- 4 ) 



¥ 

g 

fl 

^ g 


The value of ( 

constants C, , G , 
1 g 

n, 

m for use : 

in eqmtion A -1 

to A -4 are 

given below: 







Type of Plow 


t-t 


v-t 

t-v 

v-v 

n 


0.2 


1.0 

0.2 

1.0 

m 


0.2 


0.2 

1.0 

1.0 

h 


0.046 


16.0 

0.046 

16.0 

c 

g 


0.046 


0.046 

16.0 

16.0 


The follovdng eq.uatious are obtained by l.east square techniques for the 

^ 

parameter ip 1 I /L±I. ) I and R (fi^action of tube filled 

^ L ( ^ l/tp /1 a L/gj ® 

with the vapor) as a function of corresponding X, using curves given by 

Lodchart & Martinelli for and E . 

g g 

= Exp ( 1 .54625+0.48518 log + 0.044 (log X^^)^ ) (A~ 5 ) 

= Exp (1.56597+0.49115 log X^ + 0.05209 (log X^)^ ) (1-6) 
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■p = Exp (1 •55883+0.521 53 log + O.05963 (log 

- 0.0015 (log X^^)^ ~ 0.0004 (log ) (A-7) 

= Exp (1 .0705+0.5520 log X^^ + 0.0881 (log X 

V V V Afy^ 

-0.0026 (log X^)^ ~ 0.0012 (log X^)^ ) (A-8) 

R = Exp (-0.2255-0 *10046 log X - 0.0404 (log X)^ -O.OO66 (log X)^) (A-9) 

O 


J> ( / '^ P ) / / ^ P \ I 

A j P ^ I I and (fraction of tube filled with 


/■ pj 

the liquor) can he given hy Eq.A-10 andA-11. 


4 >^ = / X (A-10) 

= 1-Rg (A-11) 


Rirameter (flow type modulus for liquid) and /3 (flow type modulus for 
•vapor) are given hy eqtation A-12 andA-l5» 





p4 ./ 5-53 T, 2 , g 

®1 ’^Itv riw (Av12) 

= R ^ c/ . 2‘ _ 2 / 2 

g fgtv -Kg 


( hydraifLic diameter of liquid flow) and D 
flow) can he given hy Eq« A-14 and A-15* 


/hydralic diameter of 
gv 
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/- 


P I 


d 


\ \ J 


/ D \ 
/ p 

\ " / 


E 


(a-14) 


(a-15) 


Total frictional pressure drop can be calctilated by Eq,. A-16 by using 
eqiiations A-1 to A-15* 


, 2-n 


-• m (2hi) 


£l 

^ & L I 


2/-^] ^l-^l "1 


rri 


o( 


(n-2) i P 


(5-n) 


frio 


D 


(5-n) p 
) ^ 1 


III \ (A-16) 

\ D, 1 


g 




For the case when E '"n 2100 and thsre is 

egp / exp f 

correction by For this range ®avis defined the parameter X 

given by Eq .A-1 7 • 

/h / 

X . 0.19 ! ( Z — \ 

S ] \ / 


'.0.5 , i0-l *2 , 0-135 




Where Vm = 




\ 

W 


g/ 


^p^ / 


j (A-1 7) 


• f. 


f ! ' ’ 


-tr D 

// p 


\ ® / / 

Notations: 

(4 P/4i L) . = pressure drop per unit length due to friction,lb-f/ ft 

f ri c 

W . ¥ = liquid flow rate and ^s flow rate respectively, Ib/sec. 

1 g 

= Pia of pipe, ft. viscosity, Ib-m/sec.ft 

=3 Pensity Ib-m/cu ft. 
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APPENDIX B 

COimHISON OF fW mmmBS to ClLCUDiTE NONBOILING SECTION HE&T TEAHSPER 

CO-EFFICIENT 


For nonboiling section hra.t transfer co-efficient, two 
correlations can be used? equation given by Badger et.al.^^^ and equation 
of Martin .olli& Bcelter^^'^^ fe,dger's correlation gives overall ?fhereas 
Martinelli & Boelter's expression would give convective heat transfer 
co-efficient for inside liquid and includes the effect of ns.tural convection 
Badger's correlation- 



Martinelli & Boelter's 


^ ^ \b ^ \ 

I 

coi^relation- 


0.7 


(B-1) 


Nu = 1.75 FI E 

a.m. 


GZ + 0.0722 P2/ Gr Pr \ 

xTST / 


"O4 


/ w 


(B-2) 


3.14 N-u / GZ 
where FI = 


log ( 2 + 3.14 Nu /GZ ) (2-3.14 Nu ^ /GZ^ ' ) 

Q» •HI • 3» •in • s* *111 • O' *111 • 


F2 = I-.4O8 Tf Fu /GZ 


Nu 


a •in 




V 1.85 


(B-5) 


(b-4) 


GZ 


a .m . 


Fu 


. - 1.63 1.85<_-^-<2.0 (B- 5 ) 


GZ 


a. in. 


X = nonboiling section length, ft 


Over all heat transfer 00— efficient and nonboiling sectic® length are 
calculated for process conditions of ifedger date 7 (see Table 6 .I) by 
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using eq,. B-1 and eq,. B-2| and are ccapared beloTsr; 

Table B-1 ; Compariscn of Results by B3.dger & liartinelli* s Correlations 



Overall Heat Transfer 

Honb oiling section 


Co-efficient Btu/bu? sq ft °P 

length, ft 

Bidger et.al^^^ 

90 

6.32 

(47) 

Martinelli & Boelter^ ^ 

84 

6.55 


As clear from Table B— 1 , both equs.tions lead to nearly same 
results and for computational purposes it is simpler to use ledger 
correlation than Martinelli and B-'^lter expressioni in the Latter 
expression nonboiling section length is an implicit variable which will 
require a trial & error procedure for solution. 
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Table C.1 - Ar:$j||lne Eciuation Copstatits for ifeter 


Temp^Hange 

°C 

■ AntoiJJe Equation 

Constant 

A 

B 

c 

0 - 

30 

8.184254 

1791.3 

238.1 

30 - 

40 

8.1395986 

1767.262 

256.29 

40 - 

50 

8.0886767 

1739.351 

234.10 

50 - 

6 o 

,8.0464204 

1715.429 

232.14 

6o - 

70 

8.0116295 

1695.167 

230.41 

70 - 

80 

7.9845588 

1678.948 

228.97 

80 - 

90 

7.9634288 

1665.924 

227.77 

90 - 

100 

7.9485960 

1656.390 

226.86 

100 - 

150 

7.9186968 

1636.909 

224.92 


log^Q P = A - 


B 

T +" C 


T - tonperature, °C 
P = pressure, mm of IBg 
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Table 6 .2 - An’feoine Eq.TJation Constants for Water 


Pressta?e Range Antoini Equation Constants 


mm ox Hg 

A 

B 

C 

5-30 

8.184254 

1791.3 

238.1 

30 - 35 

8.1393986 

1767.262 

236.29 

35-90 

8.0886767 

1739.351 

234.10 

90 - 150 

8.0464202 

1715.429 

232.14 

VJ1 

0 

1 

ro 

o 

8.0116295 

1695.167 

230.41 

230 - 340 

7.9845588 

1678.948 

228 .97 

340 - 520 

7.9654288 

1665.924 

227.77 

520 - 76 Q 

7.948396 

1656.59 

226.86 

760 -1500 

7.9186968 

1636.909 

224.92 


B 

T = - C 

A-log^^P 

T = temperature, °C 
P = 


pressure, mm of Hg 
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APPSTDIX D 

GCTgPT’BR PHOGRaiSS LISTHTG 

C1 Cross sectional area of the tube 

TS Steam temperature 

PP Inter na-1 dia. of tti’oe 

XL Length of the tube 

DZ Pifferential length increment for computation 

WL Liquoi’ mass flew rate 

TI Inlet temperature of the feed 

S Solid content in the feed 

P Pressure in the tube 

Wr ?apor mss flow rate 

BPT Boiling point elevation 

GPELP Frictional pressure gx'adient 

PIT Overall heat transfer co-efficient 

Q Beat flux 

PI Temperature difference for the differential section 

PPT Pressure drop in the differential section 

QT1 Heat transfered in the differentia.1 section 

QT Total hea,t transfei‘ed 

X Yapor quality: 

Y Homogeneous phase velocity 

PEG Pen si ty of vapor 
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MUG 

IU7 

CP 

PX 

PS 

EGI 

EGO 

PV 

HTB 

US2 

PEL 

CL 

UN 

Xf-TUL 

WT 

EEG 

El£ 

EL 

KL 

ST 

H 

HS 

U 


Viscosity of vapor 

Lat&nt heat of vaporization for water 
Specific heat of vapor 

Vapor quality difference for the section 
Solid content change in the differential section 
Void fraction at the inlet of the section 
Void fmction at the outlet of the section 
^mogeneous phase velocity difference in the section 
Heat transfered in the boiling section 
Average boiling heat transfer co-efficient 
Pensity of liquor 
Specific heat of liquor 

Nonboiling section ovarall halt trasfer co-effioiont 

Viscosity of liquor 

Feed liquor mass flow rate 

Eeynold no for vapor 

Eeynold no for liquor 

Liquid volume fraction in the tube 

Thermal conductivity of liquor 

Surface tension of liquor 

Inside tube heat transfer co— efficient for liquor 

(Hfeill + St^im side) Heat transfer co-efficient 

Overall heat transfer co-efficient in the boiling section 
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IM Mean diameter of the tube 

DPO Outside dianeter of the tube 

SK Thermal condu.ctivity of tube material 

HO St^rn side heat transfer co-efficient 

Ant cine equation con stents 


A,B,C 




H 




Hi 


Mfn 
«.** ' 


^^i Cfr-'. .S'^fTIMFOCBf FAGI' S0*-‘ C»tHM^P ;.K.AG*r'W'L 

fiF-Cf.:M4 

FiVAPm/.ICR USIIG HP^IUGf '-inuS «CD'"L'" 
WOOD BLACK LIQUOR Iv ft TyB’*, HAVING Cl 

'common C%tS,DP 
COMMON/L; VFLI /XL 
> i: :'i COMHON/L'"VFL?/DZ 
V'vF, BEAD IftXLtOZ.DF 

''p.FftO 37,WLfTSfT l,S,P 
FOR MAT{'»Fe.'^ 3 
1;i¥;':F0FMAT(iF10.3 ) 

'/'K'iVwGFO* ' ' 

"'Pti'ftHftT t *•''•' *" *•“*•* IN Fy 7”* ' *■ 

PRINT, 'iOFt Witts* I tS,P 

'loa,' FORMAT i« ' WL«t' «*F jo.i, L8» M/S€C';iND«f /*, T! 

DiGf-::-' ,f ^'t/# ,*^olid co^nt"ni,''= 
''■'''',/» ,**FIO.?,«^ L8F/5-0 FT *> I 

'fb'tnt w 7 

' 10 ?. format { > » *r « '«' - # t. it V *,;»>»,» . sr-- m %' .- • # ■'OUTPUT—- - 




« , / «' P 


FVAPO^ATOP ^0 


'FVAPI wit wG,» 




iSE^Cft'FO' FOF BOILING POINT f;LrV.eTTfK«#‘*t* 
iT4'.l'9*S*9./C'- 

t:Pfl3:T,35t3,?. ' ^ 

^jtQFtF'GCJC 35 ,, 










DPF»«*GPPLF*DZ1 . ' " 

■'FORMAT{» PP=SSUF5f® *tF?A.8*» LS/SQ FT *t//*' TOTAL C"C?^ 

L8/SC FT ■'•-,//■'? FR-:IOTI,OWAL FR'«FSUR!5 DROP* ••*' , F: ; .8 , LS/i 
/''■FT #) ■" '' ’■ '' ■' ""' ' 

■ ' FORMAT! t IOC IF«) I '■■■'' '■ ' ' ' 

print 1.0, WL,S,T, U, UN, Q.QTT.QT .«.oco'«r 

'''..' 10 FORMAT!# MASS FLOW RAT*^ OF BLACK LICUCR= #,Flfe-4,® LB-'-M PW S,;- 
/NO «,//# SOLID CCNTFNT= *,Fi,6.4,« PFPC^Nf *,//* TEMPFRATURE* *,F' 
/,4,# DEGRgF FAHR:'NHEIT *^-,//» FLOW-VFLOCITY* *,Fia,4,« FT_^P£F_S'T 
' /0 *,//»■ OVHR-ALL Hf-AT TRANSFER CO-cFFICIFNT* *,F|&.4,'«^ «'TU/HR«'Sti 
.'/T.DF6 F *s//* HEAT FLUX= *,F16.4,* B1U/HR,SG FT *,//« H';AT TRANS 
'^r VRRg'D IN THfi ScCT ION= ''*,016-4,# BTU/HP #,//* TOTAL, HLAt TP.ANSFhKi 
"'/ UR 'TO THIS SECTION* *,FT6.4,# 8TU/HP *> ■ , ' , ' ■' 

."'N'^N-Fi , ' ' ' 

, 'Print 14, T'X , ^ 

' ''14 'FORM'AT! , #'801110 C POINT AT THIS PRiSSUR'”'* ♦,FI,;644*# D«^I1 rE ■ F/'#,i 

XL*KL«02? ' ' 

■/ i'E; FORMAT!#; RPMA INIFG HEIGHT* #,F16.2,* L-NgTH I NCR TMF'NT* Fl6'..i» J 

I 13 FORMAT!#'' COMPLTAtlON FOR 'SECTION *,1?,# HAS BEEN 'COMPLY' Tf'O *! ; 

,'''':PRIN:T, 9 ' ' ' „ * , , ' " 

' '■ , ' ' ' , , ' ' ' '' , ' " ''' '' ' 
''CONTINUE, ' ' ' ', . , , ,;' ,' ' , ,' ■' ' :' 

'3 ' ' ^ ' ' ' ' ''I ‘ \ < ' ' ' 

''';,^'\«j,j',,^,;'lp:,OpMATI». C,ALC0LA1ION FOr n'-yT L'FNGTH INCT«MpNT»’'n^('*i ‘I 

/'■■"C3Te'0'*0'', ''' ' , ' '' ' 

':■' ''G,6,;',T0 'J, '.' /, ' 

■' '■', ''''3 PRINT I. ' ' ■ ■ ' '"■' ' ' ' ■ 

S for NONBOILING •S-CtlC?! IS CCMPLST5 v'''''P^'V;, 

'/i ''"''4.^1' t '36 ' - , ' ' ’ ' ' • . . ‘ ' i> !> 

ts NO' NO'-eOlLlNG SFCTION'#!/ ■ ; ' 

?;'::0?;,,'','';'r'p'l'E^NO ','STarAICN 

‘';-KN'W^iLi^'Oif '' '' ' ,'', , ' '-,1’,',,',' 

',.:'‘'’;'X^swG'/'WT' ■' I,' ' ' ' , ' ■",' viy'' 




M'tL: :rp.« T , OT, Ce-G ,0'^ G , x nug, H V, C P ! 

|4#,0P*CZ#CTS*'' *.') /! ;Y^00--#W^*HV J 


I W t , K i .S , 1 , P ♦ V , n ' , D S , 0 T , D P P , D \ .t T P , Q , # L , «G U , , P G 





BBBh 




, ' PRINT 9 ■ ^ ' 

PRINT 28t»GI . 

f 'FORMAT!* VOID FRACTION T INLET* *,F26.'^ I , ' 

= ' PRINT 2&»WLfWG»X#T,V.,Pt:;.»UfeV®Tl,QTfCPP ' ’ 

i FORMAT!* WL= •,F16.4t» L'8-M/SECOND WO* 

' //« X* *iF|6«%,/» T® *tF.6-4-t» OEGRFf- F «'♦/* V»,»,fFl6.' 
/ *,/* PRESSURE® «tFl6.4,» L8F/SQ.FT «,/« SCLlO CCNTEN' 
/ PERCENT *t/* OVER ALL HEAT TRANSF-R CO- EFFICIENT* *f 
/R*SQFT.DFGR*=E F ♦t/* HEA.T FLUX* •,F'£.4,« BTU/HR.SQFT 
/ANSFERREO IN THU' SECTION* »,F16.4,» 8TU/HR ®*/* TOTAL 
/RRED* ETU/HR PRESSURE CROP* «»FI,6-4,® 

PRINT 29,RG0 

0 FORMAT!# VOID FRATION AT OUTLET* #,FI6.T) 

PRINT 9 

',;'N«N4l ' , . ■ 

PRINT r»t N ' ^ ^ • . ", 

'Tp!FLCJAT'('N>-XNK;%24#24 ■ ' ' - ' 

,1' ■ IFt!XN»*'FLOATf ' ' ■' ,■ 

■ 'f)I?»,CXN*FLOA»! Nt MOZ ■ ^ ' ■ , / ' ; 

'i43-IF!RG0»0-9'l '34v‘i(tUo , ' 

34 I!^!DPP*OPl» ■: ' 

' iFlRGO-i'WJP.Bf 3C, 20 '' ' ' 

kG;"',tRV3. ; ' ' , '' , ; ' ' ^ ' , ,, . 

PRINT 27 ' ' . , 

, ' ' POftMAt f # R ? S I ON,' IS AN NUi AR , 


^'S' 0 ,BRC»OTI rML,! WL f T » S t GPR LF t Ut Qf UN,D*" # 0 T 5 t DPT » DZ fYXvt,'P > i ' 
»^N|ftGYf.Krj,!i?|FNTUlf*i, AND CONTii^uITY 'ECUATICMS FOP. SINCif ; ''^YOOl 0 *#* 

I', COMMON;;"', piUTSfiDP' - ; ' -• 

U^',AU''.YRLTlWL,'T,S,GPftLF,!3PF,DZI ' ■'■ , 'i " ," 

''i':;"CALi;'r'"XD3,L!s,T,Dr-Li' • ' " ' ' z,,, ■, , 

' X'CL ! S » T , CL l - - "> i".' ;; ' ' '' ,■ "V " , r '!' ' ; ,',■ 

Z;;Zi'|j"PYW'©EL* 0 -Z'^l.o*GpftLF#oz •. ■ ■■ ■" 'mV'P" ' u:,:;: ; 

Z P' ''';'";Z''Z">"ZZ; z'- ^ 

-"XUNiTs, i,s ,uN) ' "i z Ui:/,'"; v,"r:"':z '"V'.' ' 


■ 3 . 2 . 4 *R*C 1 /(CL*-' ^ 00 . 0 *OiL*U^^CI,) 
' 8 t 7 »lC 0 . 0 *CL) 


sSiizzzz^siilftlliil 


•*0TI »-0.lUf J.#p 


M:144*C*14.7 i 



)fl/--.''l:f-i^|3-!-»4'5l4^ 


■M$M 


t'- -'4 ' Y- ■:■". r.5?‘?;iv 

; •-, ■■ ' 'v' '■'^' 


m, 


TX«tlfe5&.390/t7.S484-AL06l0^P3U^al4*9?> ' ^ ' 

i;:’;':'V P«p#i44. 0*14. 7/7€0.0 , ', , ;, 

' ,,' ptsp--0PT 

■ TX'fe9-0»TX/5-0+32-0 ■'■ '''■■;■• . ■ ^ " ' ' ", 

;' '' : ^ ■ TX«TX+4. 0 ' ^ ‘ , '' ' . '' " 

IF (TX^THt 6*6 ■■■■'■ - . . 

•'"''T ' ' IFUT*»TX»*.5) 6,6 *8 
,;, «'' DZ**DZ#.8 

'TaT«'{DTl+0Tl/2. ' ■ ' 

DT»0. 

TO S " ' ' ' 

::('b ■ ' CONTINUE , , ,' 

'X':',"', , ,H,ETyRN ' ■ ' 

'" ' . END' ' ' ' ' ' 

ili'FTC, PRlT'' ■ , ■ ' . . . ' ■ ' ^ "v '" ' 

' ''“ ''"''SUBROUtlNE PRLT? fcL*t,S,.OPRLF,'DPF,DZ} ,, ,' ^ //. 

■c;i('#*###t4#*4*#^**##«««^«*##*****»o***fr«**»**^'*^******^*.**'*'''*'*"**^’'^*^*f ***^'^ 

^i»*#**CACOUTES' NOt^BOlLIHG SECTION LIQUID FRICTIONAL' PRSSSuRF-, ORQP#A#i»'',, 

' '^O'MMON 'Ct»TS,DP ' ' '- . 

i'C^LL- MOL? $,t,XM'tL> '' 

,.CALL -XOELlStTfP'a) 

"" , ''i 'C„ALL MyLtS#TS,»?FUL) , . " ^ 

■',-'¥WwL^,'UlL«Cl. i 
v',':''', ',f,Ft*D,EC#V^*0'P/XMUL 

‘'/''/.''""""IFlRE-^afOO'^O'lXf lf'2 ' ■ ' V" ■- 
,,v',i:l; F«i4. 0/PE ' ■ ^ 

,''';'--Vp^^F,«{WM0L,/.nMULH*O.;;:.^ . ' 

3,, , ' ■' ■■ 

'"A*4'.0#AL0CT0t2-.C/0.0OO. ' ; 

r- ';, '',F'^t5.v.O/«'A+'2f?8 ) J f*2.0 
" Pr-Ft't'WMUL/XM0L>«*0-l? 

:i)?»E*3*i*2wti*F*oi«v*#i2*oi<^o,FL/{DP*32**! i 










iteuTiNt',XUN?TS,T,S,WMJ . , ' ' . 

■% %>. #####^‘^ #### t# 

i»^-8i|3,l,:Ll,',NG.-' SECTION 0¥;,ftALL, HFAT COEFFIC 

A *»t 4 * # »«■#***■»•«#»' « «»#•##« -jf * 

;’'|;'.XGLrS#r,Ct L ^ 

jiC»''',NCJLi"S,T,X'MtL) , ^ , ■ ." i'',') 

lL'::iX,0PL,tS,T*O!L> ' ,' 

XMUi ■' ■ - 

F0'^:06#~lC'i#0£Ll*OEL*#%O/XMULM*C.7 , 







' '.SUBROUTINe C!^C( Wl»X,S»TtP,VtDXf'0'SfOT*DPr',DV»T»?*0f'Si<L,-.'iG,.^i- GT , 
i" ' tOHMON CUTSf DP :■ ' 

' COP^<ON/L?'VP12/DZ 

'CAit P»F< WT,X,T, S,DX,0T,:0St6PR-F,RG»PJ 
; C'ALi XOEL? S,T,D!=U -■ ^ 

' CAtL PPtT, OT, CgG,D€Gl»XWG,HV,CPJ 
IF{X-0.)J?2,?l,22 
RGI=0. 

^GO TO 23 

CALL HOLOUPfWTtXtTf StRLitRGItP) 

'' OET®ll.•RGl>»OeL^RGI*OfG 
TO«T-i-DT 
' XOffX+OX 
'■ ,S'0«»S-t-DS ’ 

CALL HOLpUPi«T»XCfTO,SO,RLO,RGO*Pi 
:'‘'c;all xo^lc s.d,TO»cfti'J . , ' 

; ''Of Tl *OF. L?>R&0*'0e6l . 

'Df'TMttl DET-fDP'T 1 }/ 2» 

'OPP*»D€T«*Di-|.^^€PRF»DZ«'OFTM*V«DV/32'.2 
' ''A«CP+DPP»»'7'60,.0/''(W.O#14,7> 

'^'TX>a6f6*3'^/! 7w9^84*ALOGlO(AM »*224.?2 
/, TX«9.0'*tX/X. 0+32.0 . 

ft»'#CHAMGF- CARD 'FOR BOILING POIfiT gLFVITION*##* 

;''',0t#TX+8PF:*? " , ^ ' 

:'n:CALL XtHSrtfXL) ' ^ 

■i; ' BXDtTWA'V^J'C 8C0 .7.t - 7-^^'0'P»DZ > 

Cs!‘,Lti;*'‘»Xi;‘'6L»0T+HV*DX+X«CP«'DT+¥*OV/f?2-2 + ; 






L GO 10 2 . , 

.'OALL;''' UAF<TtS»DT#tS»0*RG.U,WT»X,OXI 

iGO''''T0'' , . ■ ; ■ 

iCALL :Uf^BLT;t'fe*'p*D1 tOSr DPPf Qf U. WT* X,DX t 
:feW#lTS*T«-0T/2. J 

1 0/ 8 '• { T. '* * X » * CL WDT. X. , -C p.> rjT-^ V* DV / C ^ 2 
Qf^Dia I*-' . (CDs > 8, S 5 


« T-" 1 7.1 0 z/:77'®' .•; IT i v 


('kL* MT^’DJ'I 


>Up:iWTfXC.tO,SO,RLO,R&0,PO) 
'r'DT.'OgG ,DE'GI t X«UG, HVf CP ) 
,<'S0»T0,ceLi'l 
'RGO'l'*Ogll+RSO*DFGI 


C0i)8,ef9 






'GO TO 10 

S/'. OV*f DV+OVl » /2. 

' RETURN ■ . 

'filPFTC RRF '■ . ■ ' ' " ' ' ' 

#*#*♦*##*♦*♦**•#**#»■«»*«««««*«♦»'»''§'»«•#«# jf****--'- - ■ ■ f. 

'C»**«*CACULATION OF TWC PHASE FRICTIONAL PRESSURE DROP******* 
■C****#*****»*#»*«## **«***♦*•**■♦'•*«*»»****♦«*«•*«♦«*«**** *******‘^^*♦****’*♦4 
SUBROUTINE PR F { Wl ,X,T , S,OX, OT,OS iGPRF ,PG ,P i 
COMMON C^fTStOP 
■' " 'C0MM0N/LFVFL2 /DZ 

' TsT+DTZ2,0 , , 

X»X+CX/2. 

^ ' S«S40S/2-0 ■ ' . ' ' ■ ' ' ■ 

'CALL «UL(S#Tt XMUL> , 

i CiLi, XDgirS,T,Ogl> . ' ■ ' , . 

" , C*Lt 'E'PIT* DTt0i6*0r'6l,fXMUGtHVtCP) ■' ' ' : ' 

, Wfe’*X*WT ■ : 

‘ ’ ' Wil^KT'f'WG ’ ■ ■ ■ ' ' ' ' ■ 1'' ^ '"'i' ' ' 

'''’ , PFG»4.0#WG/f3-i-*'’»0P*XMU5> : ^ ' • ' — ' ' 

. REi«4,0#WL/?3-l4*DP*XMUL> ' , ■ : 

i;',',, „ 6A)DE6/DeL ' . ' 1, ' ■ , . .' ,, ^ 

i;;-'- „ "C^XMUL/'XMUG ' . ■-• 

r ' ;iriRl*^5«-:2;,oo^. 1 3,i ,1 ' . - ' ' :: ^ 

8ooo.)3t2ia - ■ 'i," /' 

V'.'i'a ''"f'VM*«wL/D<?L't,wr»/DS'€J/Cl ■' l" 

GOi'TO'^'^/. : '> ^ 

I.’F,0€G*2OOO-) 6»4f4 . ' ' , 

!V/r 4 '' ■' ' ' ■ • ' 

s', ’r'S'i M«^S'0R'T'i'A*'*|4'i«8*C**O- ? } ^ ' V \ '"V ' ' '■ 

''I'; ,7'’' ■' FG*tf»XPIl*54|,2’3 + 0»48: 18*.?LCG(Hj40-04^'-^'*f ALOCfHI ' 

viRLfp$/H''' ■ ‘ '' /:!' ■ 

'^'Vr'C'Aii .«OLDUPI»<TiX,T,S,RL,RG,PI " 

;VTQ;.’i^ ^ ^ 

. , jlEf » It- 2^0 ? .,,i 9 1 1 C , 1 0 

CG- 0 . 04 t . 

H»SCRTC{.Ct/CG ^*REG**t”0,a)*4*B*C I 
D^AUOGIHJ 

PS*EKP(l.3639t*0.49l.l5«D+0.0=^209*0*«n ■ ' 


;C' 3;W0|.OWPI'.«T-f X ,T t'S t R L ,P G, P I 
[^l#-i 4 * 2 ^ 0 »PL** 2*0 

J 3 tl 3 

if|;4:|00'0.,l 13, 12*12 






ami 


C6»I6.0 
XN®0.2 

H*SQRT« (Ct/CG }*REL**(0-8)*A*6#C) 
0®ALOG{H» 

>G«FXPri-35e83 + 0.52i*?3*0+0,.05'R&3*D** 
PL=PG/H 

CALL H0L0UP(WT,X,T,S,RL»RG,P) 
ALPA-RL*»4.*PL*«3.33 
GO TO 14 

13 HsSQRK A*B*C> 

CL«=i6.0 




0013*0#»3-'.0004<^0«'*4, ) 


''1'':',^'" Dsal06'{h> 

P6=»eXPtl.0T05+.5520#D+.088l»O**'2-.C026*D««3' 
-"'PL^PG/H/ ' ' ' ' , 

'v' CALL :HOLOUP(WT,X#TtS,RL,R,G#Ri ^ 

, /'ALPA*RL»»2.0»PL*«t-0 ■ , 

44 ;5?»SQRT« AlPA/RLi 

.'.p*2.#l4.»WL/3-14 }#*C2--::N>*CL*XMUL4«^5«N/’CDP« 

, G*'ALPA#*LX'%^2.)»^«*('?'.-XN> ' ' ' 

, 'CppFteP«fe' ' 

,PRC*6P'RF#0Z 

'"/.V TteT-'Pt72- " '■ ■■ 


001,2404*4 1 






:';',^ueR0Wl,?^'p.' :HOLD,upi^.'T,‘;,^7S»p.LfRG,p> , . iv 7 ; - 

k*444M**##it:4,#«-» ««««•« S"#’?- ####♦**■»•«*«•»•«««•?'■ «*«**«**«#«#4'4#'44#'4'44*'**A 

w'iCAGU’fcATWN OF ^010 FR^^C"-TOM, tM THE; TIBO ' .■ ' 

''i;4;0MM9W:'C' f'oitTSf CP . " 

c;Ait:tDSLi'5»T*DT;i.» , ■ 

/'./'"'OTTiftO*' ' ■ " '' ' 


4-#0F L*''**"' 

!-*|4S6*H J»*2“0rL*X«-»r ) 
RPlT:f,'DT,'0?HG lOFG'^ ,X!^UG,HV,CPi 

2fc«A If ■ 

^i*t»^',iOT/12,*A I ■ ■ ■ 


;'CJtf€R all' HEAT TPA^ISFFR COTFFiCtFNT F,0 






COMMON Cl,TSfDP 
C0HMON/Lr.VCt^/DZ 

m+OT/2. 

S«S+DS/2- 

f>=P+DPP/?- 

CALL MUL{StT,XMUl ) 

CALL XDELl StTf DEI ) 

CALL XCHS»T,CL) 

CALL KL(SfT,XKL> 

CALL PPtTf DT» DEG tOPGl t XMUG, HVtCP ) 

.2> 

'PRL*CL*XWL*3 6C0 ./XKL 
;wG=iWT*lX+DX/2« J 
WL*W't-WG 

,R'gL«4-0*WL/tDP*'S .|.4«XMUL) 

PtG'tsli .0*VfG/ LOP*2' ,.t4*X MU6» 

'A»^OEt/0SG' 

'B*sXMUC/X«UL 

Ml».C*U..3+3^-*OP HP'RL^^O-^ftRf L«'*.23*PEG''« 

O*i0*OP/tHV#XM-liiL* "ft'OO^'O) 

E^P*OP/ST , , ^ ; 

W*h2' ' '' 

';:Co ttt'sv ’ ' '■ ^ ' V 

„,,CAii XHS«T,0T^»H,C,HS> ■' ' . 




3 , 4 * ft#*, as* 8 






p*p#0pp/? 

wiMm ^ 


'f'^tJlP-OUTWe . UA F(T , S,0T ’ D" 

'OVERALL HLAT 

c«***»********''»*^**« .* #t «'^******'^''f**'’*.*'*^‘ 
COMMON C!-,iT,S»OP ' . . 

COMMON/ L VF12/OZ 
T**T+DT/P. 

S^S+DS/?. ' 

CALL XCLlSfTfCLJ 

CALL KDHLf$.ff'*D€«. J ' 

CALL PP('t 0^t'D66*OfiG3,,KFUG,HV, CPI ■ ,, 

WG«><T*tX+C)X/a }» I - . ■ ■ ' 

ST**35./I^5m.**37.2I 

!■ ¥aM6ACDPG»Cl»R5l ■ ' ^ ' 

^ A».012*IV»SOPTCOF*OEG/C ?t*'42.2H }*»0-5 , 

a»CL*SQftT{Da«ST«32-2/D’^l 

H«»A*8*3fc00'0»^- ci^ /■% ' • 

f, ■ FtJRflATl* H FOR AAMULAP ^LOw» ^ 



IMfi 








CALL XHSl T tOT ,H, (,,HS) 
G»2 ./H+l./HS 
'U = l./G 
■T«T«-DT/2. 

S=S-OS/2. 

return 


filBFTC XHS 

subroutine XHS(T iDT,HtQ,HS} 

€♦*#•#**«****».«»*«,»«« in, 

€****»CACULATES( STEAM 'IDE + WALDH 

«« #•11'* «*•)'■'*#» 

COMMON CT,TSf DP 
CDMMON/L«^ven/XL 
" COMMON/L’"VEL2/OZ 
v ' ■ .T«T+0T/2, ' ■ 

, ■■ DX'^o.ooa::-; ■ ' 

■ 0M^?»D- |?4? 

' ' DPO«0 p 182!O 

■ SK»26.0 ' , ■ 

; ' , R*I-/H+0X»DP/ fSK^OM.T+-002& 

' .fW»T+Q;/R . 

■ ; /4p 

"KMUW,.»-O0,0l4 '- 

‘ DfiW'WO.^ar.'^'ia* f-t * 

': ^^,;CAit;V'PP^T5,;07•^D^^C•,K8i,,4MUG 

i:;ho»,,ho*'B 40 o^o ' ■ ' 




'WS^IU/A ■; 
t«f«:'D'T/a 
RETIMn " 

'fHW'.'"' : ' 




£IBFTC PP 


PPlT»ET»'Dr,= fi,0"C‘ iXHUGtHVfCP) ' ' 

f #######<!*■ ######## H 

I ^PfSPmtWiS OF viATrR r/:/' 

# 4.# AX ft # ^ ^ 

:ff$2,.04/9,0#273*0 ' ' ' \ , ■' ^ 

»f;.0,d0#S#t-t.252::. 05«'Tt.-#?..-,P'ST7' C9«'*’**3. )/3i 

-40 '"'‘2-'"' ■/' • ■ . ■ ^ 

'i>^''4'C*35,HTj*.47fP-0T ■ 

'9?i).''4l22l'fe22'’'400. ■■ .'gg” *T»*2'-*'^*< 

'mmm- , ■ ' ■ ^ 

'.^•l 42 2f 1'6»“ T fi - 3 0 9 6 ;• 04# T + . 2 063 C'2 « 2 • • • 00l4'2* 

3B34T»*4« , 

33*'f924/I CE4.?I • ■ , 

'•/S*043 2.^C ■ 4 '2:''.:': 

t'tf.+-DTj«-7 2,0j/9.0 ■ , , '2-,r. 








P^^Mh^^wH^H 


SUBROUTINE KL«S»1»XKL) ^ , , 

c,' '#*»*PIHE-WOOD SLCCH LIQUOI*' THERWAL CONCUCT! VlTY4#«-« 
t' ' BAMBOO • ' ' 

«« 4«»*»**«’****'******** ****** **^'****®' 
T's5.0*(T«32*0>/9.0 ' ^ ^ 

XKL=0. P335426S"’ • (03425‘»S+ ( • 0011 3 55 + . CCCC0209*S,1 *1 

XKL-XKL/U49 
^ T«9.0*T/5. 0+32.0 

, RETURN 
■ END 
''IfiFTC MUL 

S.UBROUTIN'E MUL( SfTfXMULl 
,;;C<i**#***«*R*f ••»-l •♦*#W*^*****»'»* ♦•*'*** *• 

'C '' ' *4*RPJNE**WOOD BLACK HQUOT VISCOSITY**** 

. ,0'fiA'*-*‘#‘*'*******^'**'»-*'«** +******'"^******'*^****”*' 

, ’ ', T«5.0*tT-32.Gi/9.0 ' 

': A,»'0.471T‘-0.02472*T+G-?0''‘9*10.0**{-«.C2 

"' ' ' ' ,B»0'.ofe92B»*'0.'S455*'' 0.0*41-3. .01*T+0-lfi-' 

' ' , 0^0-002046 +0.3'le 2*1 0. 0*--M-4. 01 *T«0-916 

'''';'\i. "p=pO.’Stf3,f.4‘0.:0«*l'' 4-0ll’'0.fi29*10- C**C-«6 

'>'''' ' ' h/'^ I ' ' 

■', ''''' ' ' / X!|4ULa^FX'PI A+^«S+C4S*» 1 2-0l+OfrS**43*0 1 J 
:;''l'.XWL*KHOL*O.Ot06T2' - ' 

v'';VV,:'’:,/; '" t,«9.0*T/5-0+32'- C 








ROUt-INf'';XDa<S»rrOSLJ ■ 

'BLACK' tlOUOR ^DFMSITV**** 

■« +, « » # 4^ ■«■' #r. ->*■» f' 

4:0*TT^>4wC>/SffO, ■ . 

.:»|;00t '^vO'^O .4 9 * *T + A . 0 !' S 
»Df.L*'&-C6243' ' ■ 


41*} 


^INE' 'XCLC S.TjCLJ 

■wfioo 8 lAch Liftuap 


iOg4#'!}*4l87,0 
3 9 





^^^.'ri ■•;■■•. T?‘ ?. ‘1 












•1*’ *' v’. 




vL, ...■. ., i,*«t‘»'#fc-Vf «>*» *^<''.'^-''- 

6; f;i»'i*#s I fe L F; !■: V ? FH /, T n ’ u ' ' 
|J#:t»*»lAG<5ASE dlAti Llyji'^ 

'!»■■* «•!-.■ ».• fi'/' • •' 

'SmS/U-.V^lI //.L 

IJnS^cVA Pfis ^THF' --<',S’ 

iy:j,jtAO 5',frtJa>->-£*oF 
''"fjfAD r;'tWt,TS,'i 

I, Cl'i >'i .'f^^RM'^^IT'l i* '»"' ^ *** ^ 

' ' ^ '' ''p'ft, ^^iT lO'^'f i^Lf T$« 1 1 f S«' 

'i': ^ '’' ^/ ;!■ f ^ ^ F to ., '4 f ■ '. ^ r. ^ -•: ^ ^ 






■; ^ 


^,FlO.3 ,*'' 0PCiR&E„F; 

/nmm 


flu T PUT -f*’ 


FVftPOR^^T'OF NO 


tPl'WLf«<S-»5 '»TfP»NfJ 





^3 #4 

V ^^’ ■*v'* ^ 

^vSSiil.'iisui'.J^liiyi 
. ;V\>',1vi^i;a'*fe!ik4OTf 


<£^Ck^! 




^^8 


''k 3' ' N 0 '‘^!8 ni L f >16' SrZ ft « f ' 











i600.»«t»HVl 


:0T»O 


0Sf'Mi'0PP,3V,TR,Q,MLf WG,Uf P,Sl,rS3'r 


fSl+DT^-''^ 

XaK+DX 

PaP+DPP 


CALL.'XDPtCSTl'fDFL) 

;>¥*lfl,/ f’Of 1*6" » i '■ . ■• } t 

FR.xiri.i 

‘2 6|i Wt. f w G‘,i T' f T • '/'# P ».*■' » lif 'St 'a? * A ^ * DPP. 

''pOTHAff# WL« t* ^ 

Xii* ' #i|, FI t« !'-'t F:' '..4,* OfGRs 

:/■#',/#: pA:';S§*,J3^t,» LSF/Se-FT' 

/ P?;RC£;i»T,.-*f /,^^ (■XL MR RT TP^MR-F*’ 

/a, S^Pt-O«0P:u^/, = ^ ^ FLlIZ* JtF 


# |f|G® # f F'l %f #* L S'W/ S|.iCyMO^^'f' 

-,/» V* , FT/ SECOND: 

/'• SOLID CONI t' ,*tF4;6*.Af * 
y- F.FFICIFMT*. *',Fla,<^fF‘';BTL3/H 
*•«» 8TU/H'R-S'eFf'.*i/,»''M'F^r, fR 
l}/m /«- ‘.TOTttL Hf;AT 
RTipa \LBF/'SfeFT'.»^r 


k": OUTLr;l 


llit'J NM-* 












C -t l. ■ !; L ! ■ ; ',DrL ) ■ 

^ C iL U! » '^.T - fCL ) ■','' ■■■■ 

i.'= wl,/{[> - .^■- ) 

r,:..' ;■=«:; ■" ■'■■■.:■■'" - - 



C-UL ;.;u < ■=.% r,s»UMi 
c =j.-..H.n.:: T 
:-^ = r-Vl:.U 

R*^DZ/(CL'> 600.0»D..L*U*C1) 

■ F=DZ/(‘ . <'. 00.?}“-Ci» 

^ ' '■ " "' " " ' ' ■ " 

TF{ -BMD' -■■'■ h'O.i' r*»?. t 

t=t+idt: H')/ .0 

dth;t: ' , ' , , ■ ; . ■ 

■r=7+(nT /fo.o ' , ,■■■■'■,'■ 

r = r'- + D?T 

? = ?--■■• O.CV ( ii- « 'O'f: 

T ;•> n <• :■■ - . • 9.. / i . 9^ Si-' tDG; '0 ( P J 9 } -ZZ 4 . 9 
P = P>-DPT 

■^''=9.0''-" / . '-hSv ,0 ; 

TF{(T. r 'i . , . . ^ 

t:!Z=DZ-.F- ■ ' ' ' ■ •■ •■ ' " : Fv -S, 

T = T' (bT 

t- DT=0. ■' 

COMTfMu'^ 

£!3FTC P^LT • ' ■ ' v ^ 

SUBkOUTl'iF , WLfftFfbPKLFf.OPF'fOZJ ' 

^ '?'■»'' »•:•*'; f' 5 » »»*•*?- 

C»*»'#*CACUU"f;': NnMe'OUI'43 ■ S,.-T.TIO'b tWID fRIClIOFIAL P.RcSSJR? ■ 

■*••? » # r'- *• 4 f<,' H 1 " *'•■>■ -•s.r.jtj' ?4 ; '■ »■’■ ^ ®' ^*^*‘^^* *■ **' *^^' ^ ^ 

COMMON c: , TF. OP , , ■ '• ' 













'itiv 


IMi^ 


..V;,^-(-st-V!:,.:-f;.; 




SMm 




A,;.; 

g.|(gt$i;''^- 


pi%VAi#f^jfmsOllI^!G !;.Cwr:.'yHi -JV P.'LI H- ^r.-4'-,F-"R CI>;.-,FFIClfi’4T«^'-»* ' , 

' CALt'''MUtC sit S •//■'UJL ) ' ' ' ' ' ' ' , ' " :■ ’ " 

:':V'':XALL 40FH St ft -J.-'L) ^V"■.^ 

/ ' ' ■ ' '"XWJL»5600.0''^XMUL , ' '" 

'',';,y ^ ; '*^uSo!'LmCL?‘' 0.?-LT-;..0'L-^- .0/''r-lUL)'---0-? 

:;;|f"; ' "' '' ' ' ' ' " ''''• ' 

Sti 6RyuTl 1^^^^ .”^C(kri » .' t F . 7 * ■' 5 V ; ’.}:jStr>i t3Pi fDVtT^nUtWLf 

.’-'TIJM :-'!0 COr-T* HUI ly ' j*a04T I 

; ';''' %p , ■^!A?L^P^ F» wft'N rf s ,u , oi r os, -^pr f * ■■ g, F ) 

' ", ":/;V t’.V.MGtHVtCP J 






''»'j'Tf:V|TTS»»t?fSGI,?> 
i-F'Di-L +'*(Gl . : 




;;'f«;t#).rif Tit'SCi^F Lo, mh fj 


, , , ' ' ' . , 

mi/iOll-Hl *f?.GI } > + X»‘»:'2/tD£-5*HSI ).. ,■ ■;: 
ii :» «• 2 / ? OS t : 1, . - ^ GOJ J +X ?. / 1 Pi ,S3. *5^ s PF 

‘''l^,^*"'Gi^RFwOZ- rfT-' 9- : ' 

otto-*a/,tl«^‘vO.'«’'',^'-*Tl, ' ^' '■•■ ■ ■■■^■ 

iTf / 1 T- 94 8 # «- ft L CF., ' 0 UO H 27 ; 4 . % 


P&IMT '{LFy'A'TI'-.'H 









nH" -i "I 








"CALL UHlMTf *:, P,DT,D$.DP"?,Q»U,wT» 

,«*:t|'#ffS*l-t»T/ 2 - } , , ' /.V: 

'iW:i.**CQ/B“{ i ♦•>'J'--CL*‘OT'- :r-CPi-DT-A^ 


•H T# *Z / ( 2. * 3 2 . 2 ft ?7 1 - 1 7 *:& l*'»'2>; 


tF'fABS( O'*;- [>:<■,: i-.ooo’*»'6,-, , 
'FbRMAtC 60 (lH^U ) 

Ik^Iox'+oxuV?. . 


' ; , : ,wL^^wt»wG 

’*' ' ' ' ' ' ' I :MMdN'":G '»'T$,f ciP' 






f yMP»: 2 -yjDR-*S 2 *^i >*> 3 - 














f''/ ":'' D»'ALC}S(H> 

fV;) Tf^;; p^.ii^,XP 11.5 7 -t-V * ^9 

’ ALPA^RL^'j-a-O^ PLi- 
I'lb '’irreg*Jqoo.j LUr:-’r 

m-: ,tP?RSL*g000.r;7f fUf 
|'2 CL«0.0^/.- 


1+0. 0* viOfiOi-- ■ » 




*Al.06iH> 


;•■' .;+0« Q‘'''96;' 


























"pwP+DPP /'.-■’* ‘ 
'iC^Li-'WLtS^^T, KHUi > 

'tALi^ Xf)El<"5,r ,0: L) 


;aCi^X0UL#t^iOC-./;?KL' 

*WT#(X'#DI7? .'} 


w;??7t3'P*'- .t MI.SC-) 








mm 




V,.? > ' '^*' ' ’” f| 

Kv ^v-'V' '• ' vS/.!^Vr^ 4 C'-:^ 






R a ^ i v» «Q^. T { DF^-' I'f :> / r : i . :: > > j •>-- ?•■■ o . - ^ ^ ' 

SORT( £)':L--''-Sr-^-?. ■:.'VD-^'J ' 

# 3600 . r - ' , . ''■ , 

T(# H F.iK A*!MUL?R FLaw= »,F, 

';;;;dVF: ■.'■CALL, XHSITfOT »h,VTHS) : ," , - ■ 

' G** 1 . /H+ A « / HS ■' 

#;;;:,, :; 0 "i,-/ti , , ,,;... 

'•. .,,,T*T^DT/, 2 ., 

fi^FTc^xH| 

C # * # #’4 # # ## ^ * * # » * # * ' > +> ■'‘- «' ■‘‘’- > 5' ■^' f- ’■ ”• »■’ ■“* !'■ :*■ 'r 'i- S'- ■-■'»•!■?> S‘ * fr #»■####»*»##•##### 4 i 

C>*#44C'ACUiAT,nS't St3?*;H’ %iO"*Hf-LL}i''y-T T-^Kf'r;r.:F COrFFlCl '"';V'"',, "' 

c #**»##*♦#♦«#«■ r- ###,' ,H, ^ ^;q:, .*. f. .> ■:■ <■•■ *■ ?•■»• i;- %• it ■K- » w ■■■ s- •» r- <v •r* # i * #«'- # 1 - * ® #**■#*#»'##«■»###*#* ♦^ 

/' ^ / ' CO'MMpf^/LiAv'iiiV/XI \ ' ' , ' ' ' ^ , ^;',: 


i0fO'ltF;MUGfHV',GP:l' -■ . 
i;?,’?i4#2.»32.T;:#Hy7f XL #X«UW*;C 






HV*HV#45::.*=‘9?. ^r/3 054.!’. -' ' '' — ■ ■' 

' ■ .;\:;.T*9^0»T/5.0 + -5i.0 

T|,=5.0*UT+nT )-^2.0}/9.0 . 

i?(Cv '• 'P'EC/l»2l493.fl E2l'6B2-4C}0.3 567283?#fri't'2.599l6M*Tl»#2 

" / ■ ■ ■ ' 

.OEG|.=62.43/Df 
- RETURN 
FHD . 

iXBFTC MUL 

.SUBROUTINE MJL<S»T,XMUL> 

C*###*#* •Ji***##*-^ f-lf-frKi*#-*.*! »##*#«.» #*«»*' 

'Cv;;:; #»**BASGASk 8L!|CK LIQUTR SCOS IT Y#«*« 


0 >- 79 . o^i 


A«18. 276186-. 50763'9»r+.0034?i^T»»-2 " 
I *.427 42 3 +. 043 05.4»T - , 0003 074* T *#2 
C«*Ol28T?-. 00 0698*1+. 0000064*7**2 
XMUi»6XP( A+B» S'+C*S**? ) 
KMUL«XMUL*.000672 


RETURN 


crBFTC'-'KDlt 


LIQJ3R PiNSlTT#-*##-!: 


'■pl l»l 4d20 1 ■♦■ . O'Oti *S +'f i»T' 


Q6mf4^ ;.'toP6oi2*s~r 


, V‘‘‘ A 




SUBROUTlNf''Xi:L(! 




' 0 ; : 


s V \ ; 

r. 't ‘■Ji . ‘ ’ 


'•' i'' 



















.#•¥#«£ 


iJ'''fIBJnB FAGFS0?5,NAMP R-P.SHLKLA 

£IBFTC MAIN ■'■ '. 

USES SLIP MCDFL FQP 8ADG6R,S 0A7A SIMULATIOA 

r'^ BAOf,eR,S OA^a Slf-ULATinN 

COMMON C"tTS»DP 
COMMON/L' VFLI/)!L 
COMMON/L" V?-'L?./DZ 
DO E I=xtJ‘- 
P’=6D I6,XL,D2,DF 
FORMAT{'’FP.4-} 

PPAD 17, WLtTSfT l,S,P 
17 FORMAK SFIO.A) 

' PRINT 7,1 

, 3 FORMAK# BADGER DATA MO - 

W G * 0 <• 

PRINT 101 . : 

101 FORMAT! . 

PRINT lOJ, WL,TS,1I,S,P 

103 FORMAT!* WL» *,F1C3,» II 
f T»*,F10.2,» CEGFEF F *,S 
/“ *iF10.r>,» LOF/SQ FT *1 
PRINT 107 

102 FORMAT!* 


-' •< 5 c c c .. , 
^XiCCC'-C 
r-KACCC '. C 

'K^ccc-.-c 
AKA00C5C^ 
AKAGCC6C 
AKACCCTC 
AKACCQ8C 
AKACCCgC 
AKACOICC 
AKAGOll'C 
AKAOOTIO 
AKAOO,130 
'' ■ ■ ■ AKAOOIA'C 

«»*...«»# j AKAC0150 

■ AKAOCX''feC 
■06GRKE F' *,/*AKAO'0t7C 
percent ■«,/# PAKA0OI80 

AKAOO'JgO 
AKACCICC 
'*•» AKAOCaiC 
■AKA0022C 
A1CAC023C 
,.AKAqQ ,240 
''' -AKA<!0?5C 

E > vA«A 0 ,t 2 fO 
'■ :\AKACC?70 

■ ||kJI00?8C 
3KA0029G 
.■AKA003CG 

A«ACC'3Ci 
ARA0031C 
'AKACC320 
AKACC33G 
AKA003AC 
:AKAO:0^'5C 

;akaoc 36 o 
'AKA0037C 
■v AKA00.7 8G 


OUTPUT 


^Ct«II.l4*0P**2- )/4. 

.T«TI ' 

ALL EVAP! WL, feG, «,T,P,Nf3 
'ONTINUg ' /• - . ■ 




'tCJfVAP ' ' •' . 

'':;:'SII8R0UTlNg evAPUt,w6,St;T,R,;Nfl ' 
;#RSmGii EVAPORATOR tACUtAtI6N'|#A*»#: 
COMMON ci,Ts, OP , 

XOMMON/LFVPLI/ML 
,COMMON/L€v6i2/D2 . ^ ^ Ey?'/ 

’''XN^KL/OZ » . • ',;E ■/’. 

OTpo.o '■ 


TTlPtePT) 

T+.19«'S*9 


A 0 C 4 ;;C 



a pORMATC#' calculation FOR NON8QILIN6 SECTION IS CCMPlETE' *> 
R1NT.B6 ■ ■ , 

'■pORMATf* THERE IS NO NONBOILING SECTICN 
; SOILING SECTICN, ^O-l ' ' , ’ ...¥ .■ 


' 4 i‘i|ICA 0 C 7 ' 6 C 

'^'kAfiCieo 
.AKAOCfSC 
44 ■ AKA0O5C0 
:;■■ ,AKACjD$lC 
' ■■'AkACC 8 ' 2 C 
*'■ AKAOOeSC 
AKA0084C 
' ,A')<ACC 850 
4 4 , 4'IKA005SC 
kif' ' 'AKA0C87G 
;4"' AiL'Acoegc 
,4 ,,AKAc,Cf9C 
l,fj;.;,fAKA00*5CC 
^I'i/^AKACCSIC 

•>44 A'.KA'0GS2G 
.■’:; 4 ,:..;;'AKAOOS 7 .C 
44;;4-,,4AACCS<'C 
ij4;4:;->MAGCS5C 
,',' AKAGOc^C 
A-KAO'CSTC 

(■■KklCkjC 

B ^accssc 

.I'WAGICOG 






I"-' ; 


V*WT/'!'cX* 6 a. 43 » 


^ 1 *- . 

. " -Mh 

Sliiii 


OpifAAo-o 

nT=c. 

call PP ( T» DT, DEG *DEGI 
OX 00 . #?■. 14» C P» C Z* ? T T I / ( 




01 =u . 
DS«0.'-' 
CALL CM 
T=Y+OT 
;'= :+D/ 
7=V+DV 




' , ' ".... "S-Sl^iSs.. .. ■ * 4 i.r r.. j.’.' .. '"/••MA.J.o:.- ^ 




illMi 









,WG,X,T,V,P,S,U,Q,QTl,QT,CPP A 

ID FR ACTION AT INLET= *,F16.'') ' • ' / 

- *,F’6.4,* lB”*M/SgCCND «,/• hG= #fF16.4,# LBH/SttCOND*'* 
. 4 ,/* T= 6.4,* DFGR’ir F «,/* V= *,F16.4,* FT/SPCONDA^ 

URF= #,FiS,4,* LBF/SQ.FT SOLI D' C CMENT = •»,F16.4,* A: 

/* OVER ALL HEAT TRANSF'^R CO- rFFICIENT= *,F|.6.4,* 8TU/HA 
SF F •«,/* HEaT FLUK= ♦fFl£.4,« 8TU/HR.SCFT *,/« HEAT TRA 
N TH^^ SFCTino*^ #,F16.4,* ETU/FR *,/• TOTAL H'“AT TRAhISFEA 
.4,» ETU/HR #,/* PRESSURE ORCP* *,F16.4,* L8F/SQFT *i A 
0 ■ ' ■■ ' , ,A' 

ID FRATJON AT OUTLET= *,F16.0) A: 


/ PE 
/R*S 
/ANS 
/RRE 


-KN}2I,24,24 
T( m M.)3 2, 
TtNJMDZ 
34,3C,30 

22 , 23,23 




lu. ... Uk.. .. .1 .Jk.. JIK .1. .1..... ..Mk,, 



















-1- V.. s. !■*’ ,1 . 




-;%:V|c^‘ 




*RLT« R,T,S,CiPRLF,DPF,DZ> 

3NB0IIING SECTION LIQUID FRICTIONAL PPF«^<;t)RP DROP**#* 


^ . ,k 
liiilllslliiliii 













li?'®*#T«3600./a*14*DP#02) ' ■ ^ ^ ' ' :• 

■»|l>'C4|.t ."•?? >#CL*DT+H\#DX+X*CP»DT+&g*WT#42./ r2-*13p2*TT®-i' 

-jv^ ' J'V- ' ‘ ■ , . ' , 

ro » I rr ^ri‘ 9 '’ ■ " .•'•>''#¥1 

'M§;CAI-V'I''^AP*'T, S, Cf * i:S,efRG»M»®trJE»DX-i 

ini* I'd' ‘3 ' ' 'i ‘ ' 

¥||c 4U :uKBtT*'S, P'i'Dlf OSt0PP»Qi'O*1*t»X,C>)J,J- ■ 

f 4 ^ ox'i*'< d /$** i .'li-j 14 c t 40 t»“ x«df w:' 4 «|;.f if - *32 wf* 1? s . 


AKAO? 

AKAC/ 


l-CZ/'^AKAC.’^SC 
AKA02^^CC 
AKA €?.'%: C 
AKA02-"‘:'C 
AKA024:5C 
AKAOa^fAO 
AKA02M5t 
. 'AKAO^^ec 
■ AKA0247'C 
'AKA 024 PC 
■I,:;.;': , ,;AKAfi24SC, 
-'AkACS^CC 
>::;;;:/■ AkAC?5K^ 
ij:;',; .;vAKAca52c; 
.g.':''";. i: ^AKA0272'tJ 










COMMON C 

COMMON/L"VEL?/DZ 
CALL PPF{ WT,X,I, 5,DX,nT,DS,GPRF,RG,P} 

CALL XDFLI S,T,DEl) 

CALL PP{T,CT,CeG,DFGl,XMUG,HV, CPI . 

IF(X*’0-ir:7., 21,22 , • 

RGI=0. 

GO TO ?3 

CALL HOLOUPtWTfX iT rS ,RLI ,PGI , P) 
nHT = «l.“4GI >»ceL4RGI»D*=G 
TO=T+DT 

xo=x+ox 

so»s+os 

CALL HOLDUPf WT,XC,TO,SO,PLO,PGn, P) 

CALL XDELI S0,T0,CHL1I 
0ETl=(l. ^RGa)#DEll+RGO«OEGl 
DftM*(0ET+DFTll/2. 

ACIatl.«»X|#»2/IDfL#n. .-RGI 11 +X**2/ I DFG»RGI I 
AC0s(l.-X0|«»2/{ CELl* ll.-RGOI I +XC**2 / fOFGlwRGO I 
AM«ACO*ACI 

OPP»“»OETM*DZ-GPRF»DZ»WT»**2*AM/fCl*«2> 

A*P+DPP 

CALL VPT(A,TX) 

8P€»p19*S»9*/ 5. 

' DT*t3<+BPE-T 
’:CAtL'"fCLlS,T, CL I 

’ ?’4i!|i»X4t3/, IC |G'*# 2 « R G I <» * 2 I + f I - - X I • *3 / 1 C EL-»« 2* a - -R 6 11 **2 1 


‘AkAC2Ctti' 

AKAOSCGC 

AKA07' CC 

AKAC2" : c 

AKAC2' 2C 

AKAO?" ?C 

AKAC2’ 40 

,AKAC2V1C 

AKAC2'" 6C 

AKA02I'’0 

AKA02? 8CJ 

AKAC2 * SC 

AKA027CC 

AKAG221C 

AKAC2'’2C 

AKA02*2C 

AKAC2'’4C 

AKAC2'’5C 

AKA02^6C 

AKA02?‘’0 

4KAC278C 

AKAC2'"SC' 

AKAC2' Cn 

• 1 c 

•.TO 

AKAO-^-nc 

AKAO?"*0 

akacz^'H;' 

AKAC2'^40 

AKAO?‘'70 

AKAo..r'ec 





'"■. .' ^ ^'■','’^~3^'-S‘ -“‘i-TsA.'- 


AKMifAC 
AK#C2€E.C 
AKA'Caeto 
AK'AC2fi7C 
AKACaf f?C 
AS AO 2'^ SC 
ASAC27'CC 
A'SA0271C 
'■A;KA0272C 

AKA0'27AC 
ASAOatSt 
.'ASA027|,o 
.AKAC217C 
^AKAC2"6C 
.AKA0-27SC 
AKA02SCC 
ASAC2'SlC 
^ASA0262C 
^.:>AKA025BC 
'v:aS:'AQ2SAC 
■AKA02a50 
.■'''‘IASA02?6C 
;AKAC2S7C 
■^. "ASA 02^80 
AS AO 2 esc 
.iAKAOZSCC 
■ 'AKA-02SIC 
■ AKAC2S2C 
fi<«02S3C 
. d.<ac2S4C 
AKACaSSC 

ASAOa^TC 

',t"AC2Sg<: 

^K'scassc 

ak-.aqbccc 

■' AKAC3CIC 
''.<AC?C2C 
AKAC3CiH 


5'''';v,rAlPA=Rt»*'4-'*P L** 2-33 

f/; to 14 .' 

i$^\'tF'CRfG-2000.U0, ?,8 

<l!,t;^iIP«RiL*2000^’>8,lC,l0 ^ 

;Lv'?'Ci»I6.0 ■ ■" •' ' ‘'-V 

I S CL /C6 J *R 6 G** lr^O.8 J;» A*B^C 1 , ' 

-i||‘^'^|if|p:ll!tfe397+0.4S(Ui«0+Q*.Of2O9»D*#2| 

CALL holdup ( WT, >; tTtifRL ,PG,P) ' 

AL P A *P. L^'S' 2 - 0» F L» « 2 - 0 :•■' '= ? 

■" ^ GO ' TO 14 ■' ‘ • " “ -■=' ; ■' 

O"' IF(ReG-2000IUitl3» 1 3 

i TF('>^ -L«:- ioo.n;?, .^2,1? ■' :•■. ,..: ,■' '■• 




iiMlw 




*1 






CT/.'.O , , ' 

:'~f+cy./-, 

S = S + DS/2'.0 ■ ■ , 

CALL MUKSfT, XP'CL > ^. ;■; : •■ ^ ^-V: : 

CALL KD1:L( S,T,C'l ) - " 

CALL PP(T, CT, C6G ,DFG1 tXHUG,H V,CP} . V" ;: -I" 

WG=;;'»’'WT . ' ' • 

WL =n-' WG • 

PFG=*.0^-WG/r-i .iA«DP<<-/MUl- > 

P r. L =A . 0#W L /(■■’. -r D P-fr MU L ) 

A=WL/WG 

8 = CFG/Di!L : ■ . 

C-XMUL/XMUG 

IF(REr> 2 ^ 00 -) 3,1 ,1 
IF!REL-finOO.) 3,2 ,2 
VM = {WL/Of'L + WG/DFLl/Cl 

M».19®A«fr.9'«*'B««. l-»-« VM»*2./ tDP*32.2 ) J««0.185 

GO TO 7 

'lF(RfiG-2000*> 6,A ,4 

IFtREL-2000-) e,5 ,F ^ 

H'«SORT(A»»1-8«8*C»*0.2 > 

PG*FXP{ 1.5462 2 + 0 .4P51 8*ALQG ( H > + 0.0444* ( ALOG {H H*«2 > 
PL«P6/H 
CL =0.046 
''XN»G.2 

Ota! 41 rt'Gf M 1 ■ 

OALL HOLDUP(WT,X,T,S,RL,RG,P) 







D=ALOG«H) 

PG=F'XPU .07C5 + .5 ^?C«•^ + .ca8l«D**2~.C026«C^^«3••'-.C0l2«D«■*4>:" 
PL=PG/H •' . , 

CALL HOLnuP{WTtX,TiS|RL,RG,P) 

ALPA=RL*«2.0» FL^«2.0 ■ ■ ' 

1* f *S(,RT{ ALPA/PL) ' 

F = 2 .*(A-.^WL/3.iA )*»t ? .-aN)*CL#XMUL««:?N/ tnP««l5.-J<A)*nEU«‘32.?.) 
G = . 6LP A*» (X M-?.-)* M) 

GPP,F = F«G 

PRD=GPRF«-OZ • ^ 

T=T. DT/2. 

X=X“DX/2, ^ 

S«S- CS/2. ■ ; v; 

RETURN ' 

END ^ 

IRFTC HOLDUP ^ 

subroutine HOLDUFI WT,X,T,S,RLtRG,P) 

«#*##Caculation of void fraction in 'THE tlbf««««« ^ 

COHMON C!,CZ,TS,CP ^ 

call XDEL( S,T,DEL) 

01*0- 

CALL PPtT, DT, deg fD=Gl tXFUGfHVtCPJ : : V 

'' ;Aa«“I-«(0’-G*{| .-X J««2-0FL»'X**2I , ^ , 

K H®wOlG‘*fl.-X)**2-4.*DEL»X**2 ■ ■ '' 

r^Cfl-'tOELoX^tZ ' ■ 

RG> i«t«w0«.O) /f2-*A J ' ^ 

; IF,«tG-l.S4f 5, 5 , ■_ ■' 


AKAO^;: sc 
AKA037C'C 
AKAC32iC 
AKAC32K 
AKAG32aC 

MA0I25C 
.-AKAOpfrC 

AK'A0'331I 
AKACK'BC 
AKAG32S0 
' AKA03'-3CC 
AKAC331C 
AKA0332C^ 
AKA0333C 
.AKA0334C 
Jr. AKA03^5C 
AXA03340 
• AKAC337C 
: ;.AKAG33 8C 
'AKA03390 
■' AKAC34CC' 
AKAC341C 
AKA0342C 
^ ■ AKA0a43C 
'■ AKAC344C 
AKA034FC 
. :AKAC,34 60 
AXA0347C 
iNG«4#AKAC34PI 
7' ' . ./AKA034;£C 
' ‘ A'<aC'^‘'^SC 

AKAC" HC 
Af'AG ■■FiC 
A If. AC *2C 
'iKAO '5 '-30 
;■ aKAC^sAC 

.'kag;^-cc 
AKfic nec 
A I-; AC ■■ ‘C 
’“PC 


6-fO- 











0=C«-DP/(HV^J)1MUL«''6C0.0) " C' 

fsfi.CP/ST AK#031?C 

P = ".0 AKAC3]3C 

H?=C«.22’'>«-P‘='L«*. <'*3-s‘D-s. ■». AK4C3-'?''iC 

IF(HI- 2 AKA03T3C 

1, h=H? . AKAC?.7fcC^ 

GO TO ?-i AKAp3T7C 

2 H»H| AKA03*£C 

3 CALL XHSITf DTiH, CfHS» AkAOSJ^C 

GaJ./H+l./HS A|<^A-03SCC 

'u*i./G , AKAOatlC 

T sT'-'OT/Z. ' ARA03-,fl20 

S«=S‘ OS/L AKA035BfJ 

p*p.-0pp/^. <AKA,C3l“4t 

RgTURN AKA03^5C 

rND ' ^KAC?£6C 

filBFTC'OAf ,,AKA.03£7C 

subroutine UftF(T,SfDT,OS,QtRG,U,WT,X,DXJ , ' ^ ^ ‘ , AKA03£80 

C»**##CACULATeS OVERALL HEAT TRANSFER COEFFICIENT FOR ANNULAR , FILM FL0W*Al€AC3 8 8 : 

--'CONMON C?,TS»CP , , ^ 

Sf!:;P;' CONM'ON/LEVIIZ/DZ .-.fMACBscc 

/c,,..."?; CALL^XCLtStTtCL* "i. jl;. ■:.' ’V AKACJG3L 

,;y^‘ , '■; , CALL xoEn SiT'tOgL 1 ,. 

,CALL PRCff OT, t£G tpfel#XAU6,Hyf:4R^4«||^y2-^,t^||iv: ■, AKAC 5<35 l 

!y.V;f': (.n- t»<T-* 32 - )# 5 . / R- ‘ ' ' '■• ' AKA 03 G 7 c 

??; ST*75-83n7Voa48||l»Rfiy.C|Q02lRjrW2>>f^EM|^t:M%;^M^Iil»^^^^ •». ' AKAOB'SCC 

sT«ST/{ 4$3 ^a» 32/ '''■ ’ ^ 

-VAWG/ coE!&R.ci*,R6| 

'B»CL#S<^Rt1 BEL*lf 

' '7 ■J.f 1 ' , 1 „f 4 jsi M 8 A 3 S 0'0 » O* 4 * 

C»*«*»CACULAT£Sf STEAN '' I IO£+WAt\jH'^|;yX^^AN$fiRjcqEfFlCTFNT#*** 
r. F'lRf-'/m- H POP MfNULi’R PLQW= '*t F16.<, ) .,. . ; : ■ „ 


ilfiil 


£I8FTC 




WG^i-iT-K >: + 0 :</? . ) 

WL*WT'-WG ■ . . ... .'.■ ';■ 

RFL=4.0*WL/(DP*3 .lA^XMUL > 

FFG-* .0#'«r,/IDF*:j .lAe-XMUGJ 

A^Dfa/DEG 

B=XMUG/rMUL 

C=XKL/DP 

HK,*( 1." +''i9.-*DF )4<.PRL«*C.9#-FFL*».2T#FSG««.T4*A*», 
0 =C« DP / ( H UL^f' ' 6 C 0 . 0 ) 

f=F*'CP/ST 
F = A- 2.0 

H?=C«.22''«-PPL«*. f'G-s^D^. ■». 0*P»*.3l»F-»*.T-'=» 
iF(Hi' Hi:n,»;.t 2 
1. H=H? 

GO TO 3 

2 H»H| 

3 CALL XHSIT,DT,H, CfHS» 

GaJ./H+l./HS 

Usl./G 

TsT-DTAZ. 

S*=S‘ DS/2. 

P*p. OPP/7. 

RETURN 

END 

2IBFTC UAf 

subroutine UAF(T,SfDT,OS,QtRG,U, WT,X,DXJ 


2 5*- 8 








V :’' ''';■ 




■./.r-c 


1 

SC 

'■■ ' f c 

■ CL 


AKA 03? SC 
AKAC39CC 
AKA 0391 G 
AKA03S2C 
AKAC3S3C 
AKA03S4C 
AKAC3S5C 
AKACSSfeC 
' AKA03S7C 
: AKAOBSCG 

AKAC3SSC 
AKAC4CCC 
G AKA0AC3C 
..yd'*- AKA04C;iC 
:'''d' AKACAC3G 
AKAC 4 C 4 L 
AkA04 ? 4'. 
•; A'KAC 

ara'o 

AKAO' 

A'KAC 
AK'AC 
AKAO 
Ai<A'C 

AKA 04' ;;;c 
; !■■ c 







■■' JJjEj' r *< U ; t'' r 

'tV ' > . > 'v * • T« ; 








■ AKAC4'/4 7C 
AKAC 44 ''eq 
AKAOA^'SC 
'’AKAC-'r*CG 
AKACASIC 
/>K£C''r: c 




l>«(t 

fcT “2 

w7+.0004 5«T ^• 
ti^O' 

252ll-';p5#'f*#l 

' ’ ’ ’ ' / ' J. ,v ; 

i- » , ' \ ’ 

»,.-.85f7e-qf4 

»T#*3. )/iar|':;' ;;,:;:v, 

' 'i'.;, . : '' ',1 * 


C 4 ''»'«. C 

0^-' i c 

MUG* 

191.9+0.351*1 





0 A 1 c 

















i':' - >«3f i r - 


■ vC 

AKjSC^'^FO 

AKAC41CC 

AKAOAIfi 

AKfl0477C 

AKAC4’’8C 

AKA0479C 

AKAO^FCC 

AKAC4R1C 

‘AKACArZC 

AKAC^tJji 

AKAC^RBC 
AKA0484C 
AKA04aSC 
AKA042fcC 
AKAC487G 
AKAC4f7l 
AKA04ae0 
AKA04f 9C 
AKAC49CC 
AKA'0491G 
.fAKAC4S2C 
;^jljlAC4 43t 
:'?|KA04 ?4 C 
'■AKA04 9‘"C 
iKACAcf-c 
AKACA47G 
.5KAO‘'<;8C 


IF tP*iT,340,. AND .P.6T, 
TF(P.LT.<720.« ANO.P.GT-: 
. IF( P*LT,^760., ANO*,P^^sj*! 
lF(P,GT-760%i BO .tO''''|.d 




T = i.-*'{T“’ 2 . )/ 9 . 

XMUL = I .C 3 fc 02 '^ 4^»1 + .OU 006745 -»T ««2 

XMUL=XMUL».OOOfc7 ? 

T=9.«T/5.+72. 

PFTUPN 

FND 

^IBFTC XCFL 

5UPFOUTINE XOfU S,T,D^-Li 
C ■ WA*rr.-R 

WATirR OEMS ITY«««- •» 

0 FL = . 97 ** 2 . 4 ? 

RETURN 

-ND 

£IBFTC XCL 

SUPFOUTINF XCH 5 ,T,CL) 

C WATC-R 

C##«*»WATr:R SP CJFIC H 
CL=l.O 

return 

EMC 

£I8FTC VPT 

subroutine VPT(P,T) 

C«'»***PBFSSURe TFMPERAIURE EQUATION FOR ^ATEP«^l« 1 l 
PsP«^740./( X44.*14.7} 

IFfP.LT-30-.AND. F.GT. 5. ) GO TO 2 
1FIP.LT.35.-AND. F.GT.3C. > GO TO 3 
^ \'|F(P-LT.90-.AND. F.Gr.35. ) GO TO 4 
. lFfP.LT.l50..AN0.P.GT-9C- J GO TO 5 ■ 

'1 ‘ l‘FfP-1 iwn 0 ilT 1 Krt % - 



C = 5?;S.97 
GD TO 

5 IS, = ” . 9 & 7 *!» ? B 
B = 16to5-9?:^ 

■ C = ?."t.77 
GO TO I 

9 /s,=7. 9*8396 

8=1656.39 
0=226.86 
GO TO 1 

]0 fi =7. 91 86968 
8=3636.909 


C=224.92 

T = t 8/ < ''T-CGl Cl P n 
f r:9,«T /5. + 32. 

RETURN 

END 


fCNtRY 

20- 

0.066 


0.5 

209.'^ 


0.141 

I , 16 2 


0.0 


o'a f 

II& W V 



■' I t t 

AKA052eO 

AKAC579C 

AKACT-'CG 

AKAC5"'1C 

AKi05''2C 

AKSiCS33 0 

AK^tE34C 

AKAC5750 

AK,AC53eC 

AKAOSllC 






mri «»j « 


, ^ p T *^^ *''* *' ’ •'■'■ 


i- rntmii . *fe5 *» *» *« w ‘ I»J* Wf-i « ' 1 «« * ' * 


feL= c.cee iB‘F/srcoNC 
MjS = 20 g. 7 f CFGF 7 r F- 
^•= 16? -1C CFGR^F F 

SOLID COMf;M= O.OCC PFPC^NT 
‘:''P=123C.CCC LBF/S'C FT 


mrnm-tmtmmiaim*- m^m- * •■ 


- CLTPUT--* • 
L' NGTH H'C c 


MFNT= 0 . 2=6 


REMAINING HEIGHTS 17.7* 

CACULATICN FCP NCNRCILING t? rrvPL''-T'= 

THERE IS NC KCNBCILING S‘C’■TC^ 

AVERAGE ECILIHG H''AT TPAMf-'''® rr.’'FffGT f:T- *70.90 _ 
HEAT TRAKSFI'RRtD P-! Th*- wi'i. |f.{' .\ <"r T a = ^*9' .'‘'7.89 

VAPOR CCLLFCTFOs p8..FSUR-= -'h.Q.F 


472.41 








JEJOB CGGl':?tTTM'H10e»PAG'’'Sa?0tN4MC A.K.A6ARWAL 
‘■Cli'feTC MAIN 

= ********* 

C»»***THIS USI=S SLIP MnOCL FOR BADGRRtS DATA -SI 
Cii**»»THIS PROGRAM 1 F JSFS BDFLTFR '•■CUATION F 3 R 


COMMON C?»TS,DP 
COMMDN/L -Vri I /XL 
COMMON/LSVFL’ /DZ 
RFAD l6fXL»3Z,DP 
P.5AD 17, WL»TS,TI»5tP 

16 FORMATC “F3.A > 

17 66RMAT(5F10.3 ) 

WG-O. 

PRINT 110 
llO FORMATUHIJ 
P9INT lOl 
I 66RMATA*----** 

P9INT 03, WL# TS,TI tS,P 
1 3 66RMAT{» iilL= *,Fl().3,* 
/ T®'*fFl0.2f* DFGREE " 
/?= *fFiO»3»* _BF/SQ FT 
PRINT 102 

102 FoRMATC**-**-— T*— 

ie' '■ 

Cl*I3^l4*DP**2. J/A. 
T«TI 


INPUT 


*,F10=3,» DEGREE F *t/* 
*,F10.3,* PERCENT *,/• P 


OUTPUT 




IVA-PORATOR m. lr 


PORMATC**-—*-***"'*""*** 
it eSAPlNt,W 6 ,S*T*P,NEJ 

'STOP 




t»*#«&* 4 *&&**»*******»*f»W»»*» 7 *f*-**^lTS ' T’ 

;t»*#SlWLf EVAPORATOR CACUi->ATIOHS»Jf^^^'£^£g^^^ 

»#*»*###*****»*«**#***»R»»R**** 1 »*»#|?.^»|!'*^ttf^fff- 

SUBROUTINE ^ 

;0'|, .COMMON Cl,TS,OP’ 

.,,7 C4MM0N LEVILI /XL ’ ' 7^. 

r::vri', coMMON/LEvgtt /02 ' 













• THae IS .^3 Koi^soi Litre’s* 

SICTION.MD-I ••' '•/ :%i''w.. 


P=P+DPT 

0PF=-«G7RLF»n7 1 

PRINT eiP,DPTtDPF ' ' ' " ■' ' 

^ /•=F|fc?BlI LR/SO^Fr**//»^FRl-TTnNfP PRESSURE OROP= 

/ FT *J frictional pressure DRnP= »*Flfe*S,* lq/sq 

&0RMATii.lOUH-n 

P9INT lOf WL, S»T,u,LiM»O.QTl,QT 

•.F16.4., lb-h per seco 

/T.. WGREE //. ;^n!; ^’^nrT. TSNP5 « RT UR £. .,F15 

/4 4»//& nvFP-ALL H'l’" Tp * ^K pep r o cppteT-L PtFp6.4,* FT PER SFCON 

799NT 14, TX . ^ ' ■ • 

^ pressure^ .,FS6.4,. DEEREE F .1 

P9INT 8,XL,3Z| 

‘ length increment^ .;fi6.s.. 

= for EFSTTEN p,I3,. has been COHRLETED .) , 

IFtN^LT-KJ -GO TO 51, ' ' ' ’ 

if(n-eq«k» go To. 5? ■ . ■ ^' 

; PRINT '53 . 

FORMAT! IHI I ' ' ^ ' . . '- 

l!'*’ 'K»X42 

CONTINUE ' ■ ' ' 

lFfT"*TXI2,3,5 . . . : ■■■-■■''■■-■-■■ ■ ■ 

■; 'CONTINUE' ', ' “' " ''''' 

. ,I9INT' I ' *■ • 

i “‘-=‘'‘■*”0'' forNpxr ;^ 

,TR1NT 15 .■; 

C^i-CUiATI3M FOR ' NONBOILIN 

' f ■ ' ' ' ' \ 




31 DT= 0 , 

DS= 0 . 


X=X+DX 

V=¥+ 0 V 

p=p+Dpp * 

S=S+DS 

CALL XDEL (S,r 
V=WL/{DPL»C 1 »-{| .<>PGD) ) 
0 TI=Q* 3 . * A*Do»<i.Z 
QT=QT+QT 1 

us»us+u 

HTB=HTB+OTl 
ST=ST+TS"T+nr /?. 



PRINT 9 
print 28iRGT 

2B FORMATt# vni3 FRACTITN AT INLET 
PRINT 26t WLf WGt X f T f ¥, P t 5 » Ut 0, QT 
26 format?* WLa *»F16.4,* LB»M/SfC 
//* x~ *»Fl6-fc,/* T= »TF'6.At* 0 
/ *,/* PRFSSIRF* *,FU.4,# L8F/S 
/PERCENT *,/» dvcR all HEAT TRA 
/R. SQFT« DEGRE" F *»/» HEAT FLUX= 
/ANSFERRED IN THE SECriON= *,Fl& 
/RRED» *tFlA,ii.,# bTU/HR »,/* PRE 
PRINT 29,RG0 

»>. 29 jSORMATf# VOID FRATITN AT OUTLET 
f*., print 9 

M*N +1 • - ^ ■ 

■ R-RINT 13, M‘' ^ 

ifcnpLt.ki go TO 54 : 

' <’■ ; • IFIN.E0.KI GO TO .$$ '--A 

..■aSS .PRINT $3. ' 

■ ^4 CONTINUE w' -V: 

• . IF(FL0ATtN>»KNl2lTk'4,24‘ir'7;;j[|i.fJ^^^^ 

_ , ,2l IP ^ ^ an* FLOAT? N 1 1 •»!. J 3|if''33't$'3'’' 

: 32 0Z'*?XN-FL0AT? N> l*OZ 
; ■ 33 IF{R6a.0.9I 
■■ 34 IF(DPP-0PU 22,23,23’'’’- 

22 I F ( R GO- . , r- 0 , 0 7.7. ' 7-wC% " 

/■30 TR=3, ^ 

'■ PRINT 27. , 

'' ' ll nn * * :, REGI on' is ANNULltfPl^^ 


= LBfl/S£COHD#f 

■ v= *,Fl 6 . 4 ,# FT/SECOND 
ILIO CONTENT^^ *,F 16 . 4 ,» 
IC 1 ENT= »»F 16 , 4 ,* BTU/H 
TU/HR.SQFT *,/• HEAT TR 
»f /» TOTAL- HEAT -TRANSFE 
*tFl 6 - 4 ,* LSF/SQFT *| 




PPWMIw 

- 7 . 7,7 


->51 > 7 


fotli^'g s ten on 














C##t*#CACyLATES HO'MBOI,lIMS Sf CTfOW' OftOP*4##„ 

' ' SWBROOTINE -PS^LTf WL .;.-^'‘ 
C9“HON CT',T?,DP' ;=.^.- ^ 

CALL ,f^!UL(S,’-,KMUL) ' ,'' 

, , CALL- xD^Li s, T,Dft j -- i ; 


fIBFTC EML , . ' ' . 

t^*****'^*'**^****'*'*"**'^*^*'' *’*’'*■ *‘*'^’'*'^*^* * *** ***'**^***** **■•*****•■*■■* *'**®*'*’*'*'**'**’*»*H»»l‘' 
fyf^**«.f;|^FI?GY,MOMGNruM A»D CONTINUITY F QUA T IONS •' FOR SINGLE PHASE LIQUID#*#- 

C#****####**##*.* ft *#»#«•»#*»■*»» «#**«#»®*»»**»«.#».ii, 

SUBROUTInF F'^L!WL,T,F,&PRLF,U»C,UN,DTtOTl,DPT,OZ,TK.PJ . 

'..COHMON CUTS, DP 

5 CALL PRLT( WL, TTS,GPaLF,OPF,DZI ^ 

CALL KDFL(S, T,DtL > . ■ „ 

CALL XCLf S,r ,CL ) 

U=WL/(0El»C*O 
DPT=-DiL»n7*- 1. U#G'’-.LF«-'^7 
DELT=TS"T 

CALLXUNC WL,T, S,OT,UM3 
Q=UN»DELT 
R»=0P/?..0 

A=Q*2.0»3.I's-«R*DZ / « CL*a TOO- 0»D~L*U*C1 I 
e^oz/iT.TeiT^ioo.OffCL) 

DTl*A-8 

IFCABS(0T-0n J-O.Ul,!,. 

2 T«T+<Dn-DTlZ 2.0 
OT«Dn 
GO TO 5 

I ,T«T+IDTI+0T> /A. 0 
PwP+OPT 
CALL VPTiPfTKJ 
P»P-OPT 
IFIfX-Ti7,6»S 

t IFUT*TX)-.5j 6,6t8 , : ’:.‘5 . ^ ' ' 

0|WD2#^8 . •', .. ■■U..'. ;, ■' ’ • : : ' ^v: , /■ 

4r T»T-fDTifDT»/2. ^ : ;t ' u , -v 

• 60 to 5 ' .■- ■■■.-;■: : ■: vV^'T. 

CONTINUE' ,' ‘ / ■■:■■■•:''' ' ^ ...-' 

end' - ■/ uvu-- '■ 






^IBFTC XUM ' ' / ' , '/ ■ 

SFCTTIM OV^.-RALL HFAT TRANSFER CDEFF IC IEVT«* 

SUBROUTINE ■ XJ NC WL , T, » DT, UN ) 

COMMON Ci, T$, OF ■ ■ 

CALL XCL(S,T» CL » 

CALL MUL( SfT, xmuL > 

CALL XDEL{S,r ,D“L > 

CALL KL(S,T,XKL) 


TF=123.0^ 

Hl==250. 

TF«'5.»(TF-3?, >/4. ■ "'■ 

TW=TS-.OS»TS 
5 TW=5.»J TW-32. )/9. 

BETA»(. 000161 2»(TW.TFJ + , 000003 l»«TW**2-.TF**2H»5./i (.996 
/61 2*TF**2 + . 00 0003i»^F »»?>«•( TW-TF» *9. } 

TF=9.#TF/5.+32, 

TW»9-*TW/5.+3 2. 

10 XNU«HI*OP/XKL 

G2Mw3600.»WL»CL/ (XKL^XL » 

A=3,-14»XNU/G2M - 

B»ALOG( (l-+Ai/(2.-AM ■ • ' , ■ 

Fi*A/8 ' ' ::s- 

IFIA-LT-l. 83} GD TO ? ' ■ ''-' 

>,/.• I'FIA.tT-2. 1 SO- TO 2 . ,, ' .V*,.' ^ 

n ',' P2«o- . ' " ■ ^ '■■ , • " ' '■ '■• v=-^5i. ;■ 

7 " GO' TO 3 ' ' ' ■ ^ ■ ■ ■■ ■ - ■ ■ 

1 •-Fz-i.-.Ms** . : ■■•' , ■■;'. .. ,, :.,:::?^';:;vy '. 

- 'V', „ H2=Xf«u#xxL/» '■ ' 

Z\V' 'JP<A8StHl-rt2r;il.|.l’"GCl fo 

f.;7'?7" Hi»H2 ' . ■;'; •"'■■ 

' '"''GO .TO 10 ’ . ■ 

;:;i7'4v;-vH!PtH'l+H2}/2. '*7'\ ,, 

,:.:7:;- :r CAiLL ":XHS(T#or»H,0,-,rtS» 

'-M7 ;,'* un»T-/g ■ 7' ■• ' ;'7 .;■' 

return ;■. ^ ■ -7 , - . : :.' ■ ■ "'tf. 

:, 

FIBFTC I:.,c 





■r- 


21 RGI= 0 . 

3 , GO TO 

r 22 CALL HOLDUP(WT,XtT,Sf?L] ,P.GI,P) 
DFT=(l--PGn» DEL+PGIf DSG .. 

TO=T+DT 

?co=x+ox , - 

SO*S+DS 

CALL HOLDUP! fliT,XO,TT,SO^PLr>,«GO,P> 
CALL XDEL! SOf TO tOELI > 

D 6 Tl = {1 <'-D':-l1+?.GD«-0’'-G^ 



DETM=(DFT+nrri>/". 

ACI=n .-X)»'t;,/{DFL 4 '-r .‘•■'^GD ) +'/.»»r/(DE:G*RGI ) 

AC 0 = ( 1 .->' 0 )»-H/(DFL? i.-RGDJ |+XO»‘» 2 /(nEGl«RGOJ 



IF ( Aas< on^mi 
■' 5 • ■ Dx»i m^d%% f 
11 F0RMATC60llH«-Hi''-v^::f ' 

X -i GO TO 'iO .■ t ' 

' ' 'WG»=MT»X. ■ 

'iU';- ..t', . oS«tS#WT»C!XKJ¥L-WT» 0 X> • ,- 

: ' ; ■ RETURN ' ^ j ;%' 



£ I B FTt' 

#**#'* * 

H**' WUL AT ION P 

c#****»**»*»*»» 

SUBROUTI-^F 
COMMON c:, 

COHMON/L V'-l 
T=T+0T/2.0v 
X*X’4-0X/2- 


OP T¥U PHASE ^frictional WSSSURE DR 
F { , >' , T , S ■» DX, OT f DS , GP RF , R S t P )J . 


‘ #»»;#»*##*#•***»*** 
»«♦*#*** •.''Tv,. ^ • 


CALL rD-i(...,..-L), 


AM=ACO-+CT 

' OPP=-DFTM»DZ- GPPFi«-0?«WT •-*■?» AM/ f **2 » 

A*P+OPP 

CALL VPT(A,TX) 

BPE«-l9»S»9./5. 

0T»TX+8PE-T 

CALL XCLJSfTf CLt ' . 

AEl»X**3/t DF3*«-2»RGI»»20 + ( I . -XM'»3/ ! DEL»*2*( I - - RGn*»2 > 

AEO*XO»»? / < D= G}>«'?*RGO*»'2} + i I .-XO) »*B/ ( DELl**2*< l--RGO)'»«2> 
'AE=AEO-An 

B='MT*3600« / ( i •1 A#DP*OZ > 

C»a.-X»»CL*0T+HV*DX4'X*CP»DT + AF*WT**2/«2-»32-?.*778-l7*Cl##2i+02/T? 

/a.n ■ ' ^ , 

lFtTR-5Q.2-) GO TO 2 ' , _ ' n ' 

CALL UAF! Tt S* OT » OS f D# RG» U>¥T* X* DXT,%.^.; ; '■ .-w 1 * 4 -“^ 

GO TO 3 ' '.T-' ■■. : ^'v. , , ' 

• CAiL -UN 8 I T*S, P,OT,OS^OPf^Qf,M,f:;MT,tXj».^X| >;/ ''ys? vi;,,; : 

0 *U» I TS»T»DT/ 2* ) ' • ' ’ ■ ‘ ■'"?•■ ' ■•i r*^' ‘ - ■ 

0X1 « I a /8- jl #.f XT »C L *2 / T%:|f 2 *S? ** ^ * * *^^ "^^ 




CALL -PPfT, r>T, 0 !:(,,d=G' ,''>MUG,HV,CP r 
WG=X*WT 

WL®WT-WG 

REG«4.p»riG/ C 3 .1^»-D'P»<MUG) 

ReL*=4.0#V^L/(? .?,‘'+#OP»XMUL) 

A«WL/W7 
B=OFG/Dn 
C=XMUL/XMUG 
IF(RFG-22 00.) 

IF(RFL-8000,J ■ .?f ? 

VH=( wl/d:.l+wg/D'-g)/:'5 

H = ( VM»-B-2./ (DP*C2.1) )*»3.ie5 

GO TO 7 

IFmEG-”2000.> 6,4,4 
IF(RFL-2000-) 6,%5 

H = SORT( - 

PG=EXP« !• 5462 3+0- 485J 8»AL07 ( H> +0 ALOG ( HI H **2 ) 
PL»PG/H 
CLeO.046 
XN»0.2 
D=AL0G{HJ 

CALL H0L0UP(»<T,K,T,S,RL,RG,P) 

AL7Aa=RL*»4.*PL**3.33 
GO TO 14 

IFfR66-2000-) 10,8,8 
IFfREL-2000.> 9, 10,10 
CL«|6-0 ^ 

■ CG«0*046 
XN^l-O 

H»S,QRT( ICL/CS »*ReS**f " .0>*A#B*CJ 
0»=AL0G<8t , . ■ 

. PGw6'xra-S6®^7+0.49.H.5»-0+0-.05200*D*»2J 
PL*PG/H ■. -.-fe’',- 

iH'XALL HO*Ol»PtWT,)t|T,S,,RtfRG,f»l ■« . • , ■ 4 - ■ 

' AtRA«RL»*2 - 0* PL» *2* 0";'3v.' V/ . ' ,; 

'.H" 6 o TO 14 ., ; 

lFmEG*2000. Ill, 11,13 y' 

' iFmet:-2ooo.» «, 12,12 ;;;,4, 

, 4CL«o-04e» ^ ; ... 

' CG«16.0 

H=SCRTnCL/C''*U^s~:L»»? O.CI*A*B*Cl ■ V’ 

0 = AL0G(H) 

^ ■ P G =»6XP { 1 . 3 $ 08 3 / 0 21 0»O+'life'm^‘SM &i* 2!* -:0 j 

PL=PG/'H •• •••'■' v-Vl'* 

CALL H0LDUP{WT,X,T,8-TLf'’G, PI 


j|k ' - - 




6o 4£ *b»»'2f-4 601 2»b»' 


PL«*2.0 







P = SCRTt ALP4/^L > 

F = 3«*( .“i ■■ 

G=:^LPA«-»( XM-? . " 

GPPF=F»G 
PRO=GPR F®-DZ 
T«T-0T/2. 

X=X-0X/2. 

S = S"DS/?. 


.■,M)»CL»-XMUL**XN/ {r)P»«-( 5--XN)»DgL*'32-2 » 

/M) 


RETURN 

FND 


£ISFTC HOLDUP k«.«,«.«.*»«.#.3«.**##'»«-»*»*»*******‘*'*‘*^^'****'*'®'*'*'**'’ 

#- fc nil AT T riM nc ^nin F^fCrlON IN THc TUBfc^**** 

SUBROUTINE H 1L0UP I WT » t T* St R LtRG » P > ^ 



f, ,r ' - 


Km 




Mi 






COMMON C’,TS,DP 
CALL XDFL(S,r fQ“L) 

DT MC- 
CALL PP(T,DT, OFG.D'-Sli 
A^“2.*(DLG»(1 .”X)*»2-C 
' B=*Df:G»( 

C=2-*D+L»X**'? 

0=SQRT{B**2'"'+-'*^4»C) 

Pf ={-2+01/ (2. »+> 
R61«I-B*DJ /I2 .*A) 
IFiRG-l- »4,5, 5 
5 RGaRGl 
SO TO B 

* % ,vtF.fRG-0. IZ,3t3 
>'? i' '^RGeRGl ■ ■-: 

§ ^lli*l.*RG - '■ V'; ' ' ,: 

■RETURN ' " 

END , rv"' 

*|0RTC UN8- , ■ , ' ' / 

C*#*&#»»***»***»*^»*R®**^^*‘* 
,C*»##»iCACULATES ' OVI'RltL 'p* 

SUBROUTINE UN8(TfSiR» 
COMMON C?»T‘tO«> 
COMMOM/L~V-LE/OZ 
T=T+0T/'-:- : ' 

S=S+DS/2. 

P=P+DPP/2- 
CALL MUL(S»T, XMULI 
C 1 LL' X 4EL I S -I T ^46 L ||- '|^ 
CALL XCLtS^TtCL) , ' 
CALL KL( Sf T, <KL) ' ■ ' ' 
' CALL PP(T,DT,D*=G.D'G! 

^ f a,. 


*•;**•** 




if.’' I-' 








1 / ‘‘/I 


SftS+DS/2, -■ • ■ . ■ 

CALL XCLCStTfCLJ 

CAti X0fLIStTf46L) ■"'■"■■ 

CALL PPtTf OTf OlCtQeSlfXMUS.HV^CPI' 
«G=«T#fX+&X/?..>' ^ ”:v* '■- ,*;’ ' 

t=’IT-32.>»5./'9. '■ ,•; ■ , 

ST =75- 83076*0 - 14835S*T;--iO0O21 »T**2 

ST*ST/{4S3-Ar-?..rj 

T»9-»t/5--t-32. 

5*W7 / 1 4EG*CI y 

A'^OlZa ( Vf $0«! T ( Dp*0FS/-iSt*3E- 2'T|. I# 
8=CL»SOftT(,OFu^ST*3?.r'/opy ' . 

H8v*®#36d0.0i^^-H ' ' ;. 


REG=4.0»WG/( 3P*'5. .3.*^!5-KMUG) .. , ^ vv;?^ 

L A=DrL/4+G 

' B=XIMUG/XMUL '■,■''' ' ' 

C=XKL/OP 

• HI =C» ( i« 3 + 39. «-DP J + P5 L»*u. 9 r-R“L + * .2“' *REG«'», 34*A#«- .25»B 
^ , D=0*DP/CHV*XMuL + i£00.0» ■ " -3:3- ^ 

5=P*0P/ST :. .4 

F = A” l,« 0 

82=C Z 5'*’P9L’**.j9+D'*‘+.f9'«'r + &.*'3+F + »,3?. 

IF{Hl»H2)1,,l,P :. 

I H = H2 

GO TO ? , 

Z H=HI •» 

3 CALL XHS? TtOT 
G*1-/H+1./HS 
U=l./G 
T=T«DT/2. 

S=S»DS/2- 
' P=P-D7P/2- 

RETURN 
END 

IIBFTC UAF 

C ****♦***»****•«••»»* ■»»*»»»»»**«♦♦*»»«**#*»#*»*«»##*»#»*##*♦*#*«*#•*•#»**•• 
C****»CACULATES 0V=RAIL H'&T TRANSFER COEFFICIENT FOR ANNJLAR FILM FLOW»» 
»«•#!» *♦#**••#»*■»#**»*»»#*«•♦##****#*«■#»■»»»#**♦#»*#»*#»*»«» •t 

SOBROyTlNE mF(T,StO^»DS,Q,RG,U,WTfXtOXJ 

■-COMMON Cl,TS,DP ' ■ * , ' .i, 

^ C64M0N/LEVgL2/DZ ■ ' ■ . , --rK:'--' '"‘''’i 






■H 


■IHiRSI 


" / + . 00000" t »^T- .“-,1. ., ^....,...«„- 

''!’=HVfr45 9? 4/5 0" '■». ?? . V ,' * , - V-.' -•■ 

T=9.0* 1/5-0+12*0 "'■i't’v ' i 

Ti=5.0*(('T+0T}"»32-05/9*^ '•' ‘ ^ L**" w '•#«/■ firj,’ 4' 

45 61=21 493 - *l ,22 i 63 2*4 00.^|'‘f ||2 ^■: 

0EGl=e2.i3/0’-G1 -. '. •: ;;, r ;';>”■■■ ■■■',' '■ ' ‘".•- ■ 

■*»****♦*#**» »»»» +»+ »»#»t'#»i 4 >j|,»| i»»».£#. fy,»+.*»»»t»yf»» ♦**_ 


005 T 9602 *' 


C0!«M0f'l/L VPl'^/nz 
XL=20. 

T»T+nT/2. 

0X=-0095 
DM».I3«:5 
DPn=-X66 
SK * 21,8 - 0 

R=5 ./H + DX»0P/ {'^Kr-DM) 

P = !.-/R 

IF<Q.?8 = 0. > 50 "^U :. 

TW=T+Q/R 
70 TO 2 
TW=TS-.0''*T5 
TF*TS-3.»? TS-’TW J /4. 

7KW=. 594/3600 . 

CALL MULJ Sf T" ,XMUW) 

45W=-97*62.4S 

CAL'i PP(TS,D»- ,D = G,0=31,KMUG,HV,CP» 

H0=l- 13*( XKW» *3- »DFW»*2**32- ^.•HV / I XL*XMUW*C TS“-T W 5 5 I *»-25 
Hn*HO»3600.0 

A=DX*DP/(SK#0M)4DP/Ha«'DPn5 

HS=*l./A ■ ^ 

3=T-DT 2- 
RSTURN 
^ND 
fI66TC tP 

C;»#&|i*P8tSlCAL PROPERTIES 06 WATER VAP0R***f*,, ^ ^ 

SU8R0UTME PP«T,9T+DPS,0EGUXR»U6fMVfpPt '0 ''• ■ 

1 T*5.0»tT-12.0l/4*0+tT8.0'.* ./».■>• • ■ - ;,i0 ' I‘ ,‘u., 

■i. *T.«T.T«ooo»s.T*.45Z|po5.i«2:».w*r^ojjm^^^ 

‘xMU6=in?9+o.3»l»TT».MZE-'0T,' ■ - ■ "' 

■■ '• '' 0EGa2l493-6l?2l'63|*4l?'tt-%fipi832*t+a!-f'^^9|.63l£T|t?%«'^O05|.96024f«»3. . 

HV*=2799»i422976*’'lS‘-3‘I^'96504RT+fi206^^4|^^^**R!^’!^^4^^'^^.^T'^^*^**^"*' "'■ 







M!)L 

«• *■!»>»#■#« «• ;'r <Hf r*« 

MATi-'R VISCOSITY* 

SUBROUTT 4P MJLtS 
T=5.»(T-S^. >/ 9. 
XMUL=1.996C'’'' -.0 
XMUL = XMUL*. 000 
T==9.#T/5. + ?:'. 









B= 17 : 5 . 4;;9 
C = ?.32.i4 , 

GO TO l„ 

6 # = P. 0116 ? 9 f- 
6 = 1695.157 
C = 230.41 

GO TO I 

7 ^ = 7. 9845*^-88 
8= ! 673. 95 P 
C=228.9T ■ 

GO TO S 

8 #.=7.96347 8 6 





8=1665.97 4 
C=227.77 
60 TO 1 
9- A=7-948396 

• 8=1656-39 ■ ^ 

... C=226.86 
' GO TO I 
10 A=?. 918-6968 

8=1636.909 

ss ^ A. 431 !> 

I T=?8/U“AL0G10(P) J J 
'4'3’/ , - T«9*»T/5.+-32, . 

■V ’ ,, P»»#»144.*i4. 7/760-. 



